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INTRODUCTION 


This  report  presents  the  results  of  a combined  analytical  and  ex- 
perimental research  effort  to  investigate  the  dynamic  stability  and  trim 
characteristics  of  a Froude  scale  model  of  the  AEROCRANE  heavy  lift  vehicle. 
The  research  described  herein  is  intended  to  quantify  the  aircraft's  transfer 
functions  to  control  inputs  and  its  trim  conditions  in  hovering  and  forward 
flight  operation  in  still  air. 

An  analytical  and  experimental  investigation  of  the  hovering  dynamics 
of  a 0.1  Froude  scale  model  was  conducted  and  reported  in  Reference  1. 

During  these  investigations  it  was  determined  that  the  model  exhibited  a 
lightly  damped  retrograde  precessional  motion,  that  under  certain  conditions 
of  rotor  thrust  and  center  of  gravity  positions,  could  become  unstable. 
Although  piloted  analog  simulations  indicated  that  with  proper  motion  cues 
a pilot  could  stabilize  this  mode  with  a reasonable  level  of  effort,  it 
was  also  demonstrated  that  a remote  pilot,  with  inadequate  motion  cues, 
would  have  great  difficulty  in  controlling  the  model.  Theoretical  in- 
vestigations considering  a coupled  four-degree-of-freedom  hovering  ana- 
lytical model  showed  that  a relatively  simple  azimuthally-phased  attitude 
feedback  would  readily  stabilize  this  mode  at  very  low  levels  of  feedback 
gain.  Accordingly  it  was  recommended  in  Reference  1 that  such  a feedback 
control  system  be  incorporated  in  the  model  control  system  for  future 
experimental  investigations . 

It  was  further  recommended  in  Reference  1 that  the  analytical  model 
of  the  vehicle  be  extended  to  include  forward  flight.  It  was  considered 
that  this  extension  of  the  analytical  model  would  provide  insight  into 
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the  dynamic  behavior  of  the  vehicle  in  forward  flight  that  would  be 
valuable  in  planning  and  conducting  forward  flight  experiments  with  the 
model.  This  extended  analytical  model  would  also  provide  the  basis  for 
corroboration  of  both  hover  and  forward  flight  experimental  results  in 
quantifying  the  vehicle  transfer  functions  to  control  inputs. 

Finally,  it  was  concluded  in  Reference  1 that  certain  aspects  of  the 
model  and  model  control  system  could  be  modified  to  increase  the  safety 
of  experimental  flight  operations.  In  particular  it  was  deemed  advisable 
to  incorporate  a means  for  rapidly  achieving  a buoyant  state  at  any  flight 
condition  and  simultaneously  arresting  the  model's  rotational  motion  so  as 
to  lessen  the  probability  of  model  damage  in  the  event  of  model  control  or 
power  loss,  etc. 

The  research  efforts  reported  herein  incorporated  the  above-cited 
recommendations  and  implemented  the  conclusions  of  Reference  1 by  measuring 
the  hover  and  forward  flight  trim  and  dynamic  stability  characteristics  of 
the  AEROCRANE  heavy  lift  vehicle , using  a Froude  scale  dynamic  model  in  free- 
flight.  The  experimental  results  are  compared  with  the  results  obtained 
from  the  analytical  model. 
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ANALYTICAL  MODEL  FOR  TRIM  AMD  DYNAMIC  STABILITY 


This  section  presents  the  analytical  model  to  predict  the  forward  flight 
trim  conditions  and  the  dynamic  response  characteristics  of  the  AEROCRANE. 

The  formulation  is  complicated  by  the  fact  that  the  rotating  centerbody 
produces  lateral  force  owing  to  the  Magnus  effect*  and  thus  forward  flight 
trim  involves  a vehicle  roll  angle  as  well  as  a pitch  angle.  Further,  owing 
to  the  comparatively  large  drag  and  Magnus  force  developed  by  the  centerbody, 
the  trim  roll  and  pitch  angles  are  relatively  large  and  consequently  small 
angle  assumptions  are  not  made.  In  order  to  make  the  presentation  more 
compact  the  development  is  presented  in  matrix  notation  and  then  expanded 
to  produce  the  trim  equations.  The  dynamic  response  equations  are  obtained 
by  a perturbation  analysis  about  the  trim  condition.  The  perturbations  are 
assumed  to  be  small  angles;  however,  the  large  angle  formulation  is  retained 
for  the  trim  or  equilibrium  condition. 

I.  TRIM  ANALYSIS 

1. ) Axis  Systems 

The  following  axis  systems  are  defined: 
a.)  Gravity  Axis  Systems,  XQ  (x8,  y8,  z8) 

This  axis  system  is  oriented  such  that  z 8 is  parallel  to  the  local 
gravity  vector  and  points  downward.  x8  points  forward  in  the  direction 
of  flight  of  the  vehicle.  The  flight  velocity  of  the  vehicle  lies  in  the 
plane  of  x8,  z9  and  in  general  is  composed  of  a horizontal  velocity  UQ 
along  the  x8  axis  and  a vertical  velocity  WQ  along  the  z8  axis.  Thus  UQ 
and  WQ  describe  the  velocity  of  the  vehicle  with  respect  to  the  earth. 
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b. )  Gondola  or  Body  Axes,  XB  (x8,  yB,  zB) 

This  axis  system  is  aligned  with  the  body  or  gondola  of  the  vehicle. 

The  body  axis  xB  lies  along  the  longitudinal  axis  of  the  gondola  pointing 
forward, and  the  axis  zB  lies  along  the  shaft  or  vehicle  axis  of  rotation. 

The  orientation  of  this  axis  system  with  respect  to  the  gravity  axis  system 
is  given  by  three  rotations;  9,  the  body  pitch  angle;  0,  the  body  roll  angle; 
and  If,  the  body  yaw  angle.  These  rotations  are  performed  in  the  following 
order:  0 is  a rotation  about  the  xa  axis;  9 is  the  second  rotation  performed 
about  the  deflected  y axis  (y'  \ and  if  is  then  rotation  about  the  further 
deflected  z axis  (z").  The  cyclic  control  is  referenced  with  respect  to 
the  orientation  of  the  body  or  gondola  axis  system,  i.e.,  the  azimuth  angle 
for  cyclic  is  measured  with  respect  to  the  negative  xB  axis  and  is  positive 
in  the  direction  of  rotor  rotation. 

c. )  Rotor  Axis  System,  X„  (x„,  y„,  z„) 

This  axis  system  is  employed  in  the  derivation  of  the  rotor  forces  and 
moments.  The  z„  axisf  lies  along  the  rotor  shaft  or  axis  of  rotation  and 
the  x„  axis  lies  in  the  plane  of  the  relative  wind.  The  orientation  of  this 
axis  system  is  obtained  by  rotation  of  the  body  axis  system,  XB,  about  the 
zB  axis  by  the  sideslip  angle,  9 such  that  the  x„  axis  lies  in  the  plane 
of  the  relative  wind.  That  is,  by  definition  in  this  axis  system,  the 
velocity  component  along  the  y„  axis  is  equal  to  zero  in  equilibrium  flight. 
In  the  derivation  of  rotor  forces  and  moments,  the  azimuth  angle  is  measured 
with  respect  to  the  negative  xB  direction,  positive  in  the  direction  of 
rotor  rotation. 
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d.)  Wind  Axis  System,  Xw  (xw,  yw>  z^) 

The  centerbody  drag  and  Magnus  forces  are  defined  in  this  axis  system. 
The  xw  axis  points  in  the  direction  of  the  relative  wind.  The  orientation 
of  this  axis  system  is  obtained  by  rotation  of  the  rotor  axis  system,  Xt 
about  the  y„  axis  by  the  angle-of -attack,  a,  such  that  xw  points  into  the 
relative  wind.  That  is,  by  definition  the  resultant  flight  velocity  is 
along  xw  and  the  velocity  components  along  the  yw  and  z^  axes  are  zero. 


e.)  Shaft  Axis  System,  Xg  (xs,  yg,  zg) 

One  further  axis  system  is  employed  for  force  and  moment  resolution  and 
this  is  referred  to  as  a shaft  axis  system,  which  involves  only  rotation  of 
the  gravity  axis  system  through  the  first  two  rotations,  0,  the  vehicle  roll 
angle  and,  0,  the  vehicle  pitch  angle.  This  is  convenient  owing  to  the 
fact  that  the  two  rotations,  the  body  or  gondola  yaw  angle,  t>  and  sideslip 
angle,  0,  will  appear  as  a sum. 

These  then  are  the  five  axis  systems  employed  in  the  development.  They 


are  shown  schematically  in  Figure  1.  To  proceed  with  the  development  we 
employ  the  following  compact  notation.  The  symbol  X with  a subscript  refers 
to  a particular  axis  system  as  well  as  the  three  components  of  any  vector 
defined  with  respect  to  the  particular  axis  system,  for  example,  the  three 
velocity  components  in  that  system  or  the  forces  or  moments  with  respect  to 
that  axis  system.  Further,  the  various  rotation  matrices  are  denoted  by 
single  symbols, 
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-sin  0 

cos  0 

pitch 

cos  0 

0 

-sin  9 

0 = 

0 

1 

0 

sin  6 

0 

cos  6 

yaw 

" cos  4 

sin  if 

0 

Y = 

-sin  If 

cos  4 

0 

0 

0 

1 

sideslip 

cos  g 

sin  g 

0 

B = 

-sin  g 

cos  g 

0 

0 

0 

1 

angle- of-attack 

t 

r 

cos  o 

0 

sin  a 

A = 

0 

1 

0 

-sin  a 

0 

cos  a 

m 

(1) 


Note  also  that  since  all  of  these  matrices  correspond  to  rotations, 
the  inverse  of  any  of  these  matrices  is  equal  to  the  transpose  of  the  matrix. 
Further  it  may  be  noted  that  when  the  matrix  product  Bjr  appears  it  can  be 
expressed  as 
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COS  (0  + ♦ ) 
-sin  (0  + t ) 
0 


sin  (0  + if  ) 
cos  (0  + t ) 
1 


0 

0 

1 


thus  Bfr  is  a function  of  (0  + if)  only. 

Transformations  among  the  various  axes  are  given  by 


X,  » ♦ 0#  Xa 

X|  a B X| 

X “AX. 
w " 


xs  -e  § xa 


(2) 


(3) 


2. ) Forces 

The  various  forces  acting  on  the  vehicle  are  defined  in  various  of 
these  axis  systems  as  follows: 

a. ) Gravity  and  Buoyant  Forces 

The  weight  of  the  vehicle,  acts  in  the  positive  za  direction  at 
the  center  of  gravity  of  the  vehicle,  and  the  buoyant^  force , FB  acts  in 
the  negative  za  direction  at  the  center  of  buoyancy  of  the  vehicle.  Thus 
the  vector  of  forces  produced  by  gravity  and  buoyancy  is, 


b.)  Rotor  Forces 


The  rotor  produces  aerodynamic  forces  consisting  of  thrust,  T,  in-plane 
force  H,  and  side  force,  Y,  defined  with  respect  to  the  rotor  axis  system 
(X„).  The  thrust  is  positive  in  the  negative  z„  direction,  the  in-plane 
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force  is  positive  in  the  negative  x„  direction, and  the  side  force  is 
positive  in  the  positive  yR  direction.  These  forces  act  at  the  rotor 
hub.  Thus  the  vector  of  forces  produced  by  the  rotor  is. 


c.)  Centerbody  Forces 

The  influence  of  forward  speed  in  combination  with  centerbody  rotation 
produces  a drag  force  and  a magnus  force  which  are  defined  in  the  wind  axis 
system  (Xw).  The  drag  force,  D,  is  defined  as  positive  in  the  negative  xw 
direction  and  the  magnus  force,  FM,  is  defined  as  positive  in  the  positive 
y direction.  These  forces  are  assumed  to  act  at  the  center  of  the  center- 
body.  The  vector  of  forces  produced  by  the  centerbody  is  thus 


It  is  assumed  that  the  gondola  produces  no  aerodynamic  forces  owing 
to  its  small  size. 

The  contributions  of  these  forces  can  be  summed  to  determine  the  re- 
sultant force  acting  on  the  vehicle.  In  matrix  notation  summing  forces 
in  the  direction  of  the  shaft  axes,  the  three  force  equilibrium  equations 
are,  using  the  transformation  relationships  given  by  equation  (3), 


— ^ 

t 


>■ 

! 


[ r 

f 


j 

i 


Xg  = (ABf-)T  \ + (Bf  )T  X„  + 8*  X9  (7) 

The  first  term  represents  the  centerbody  forces,  the  second  term  the 
rotor  forces,  and  the  third  term  the  gravity  and  buoyancy  forces.  Equi- 
librium flight  is  given  by  the  condition  X = 0.  Equation  (7)  can  be 

s 

expanded  using  the  rotation  matrices  given  by  equation  (1)  and  (2)  and 
the  forces  given  by  equations  (4),  (5)  and  (6)  to  yield, 

(-W  + Fe  ) sin  9 cos  0 - H cos  (P  + if ) - D cos  (P  + if ) cos  or 
- (Y  + F„ ) sin  (p  + * ) * 0 

(W  - FB  ) sin  0 - H sin  (P  + if)  - D sin  (p  + ♦)  cos  a (8) 

+ (Y  + F„  ) cos  (p  + * ) * 0 
(w  - Fb  ) cos  0 cos  0 - T - D sin  cr  « 0 
3. ) Moments 

The  rotor  is  the  only  component  of  the  vehicle  assumed  to  produce 
direct  moments.  The  rotor  produces  a rolling  moment,  L, which  acts  about 

the  positive  x„  axis  and  is  thus  a vector  in  the  positive  x„  direction. 

# 

Similarly  the  rotor  pitching  moment  M is  represented  by  a vector  in  the 
positive  y„  direction.  Thus  the  rotor  moments  expressed  as  a vector  are 


Since  the  rotor  is  propelled  by  tip  propulsion  the  net  rotor  torque  is 
zero.  It  is  assumed  that  the  rotor  RIM  is  constant  and  consequently  the 
balance  of  yawing  moments  is  not  considered  further. 

* j 

] 

j 
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Moments  are  taken  about  the  center  of  gravity  of  the  vehicle  which  is 


assumed  to  lie  on  the  axis  of  rotation  a distance r below  the  rotor  hub.  It 

o 

is  also  assumed  as  noted  above  that  the  center  of  buoyancy  is  coincident 
with  the  rotor  hub  and  that  the  drag  and  magnus  forces  act  at  the  center  of 
the  centerbody.  Thus  taking  moments  in  the  shaft  system 

Lg  = (P*)T  L«  + HxX*  . (10) 

X*  includes  all  of  the  forces  acting  at  the  rotor  hub  and  consequently  in- 
s 

eludes  all  of  the  forces  contained  in  X with  the  exception  of  the  gravity 

s 

force.  The  vector  H is  defined  in  the  shaft  axis  system 

H = -r  k 
o s 

where  k is  a unit  vector  along  the  z axis, 
s s 

In  vector  form 


X*  = 
s 


Consequently 


H x X*  = 
s 


r Y*i 
o s s 

-rX*3. 

OS“ 


Therefore  the  summation  of  moments  using  relationships  (9)»  (10)  and  (12)  is 

EL  = L cos  (p  + *)  - M sin  (p  + *)  + r Y* 

s os  (13) 

EM  = L sin  (P  + f ) + M cos  (P  + t ) - r X* 

3 0 S 


where  from  equations  (8) 
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X*  = F*  sin  8 cos  0 - H cos  (p  + f ) - D cos  (P  + f ) cos  or 

(14) 

- (Y  + F„)  sin  (p  + p) 

Y*  = - F|  sin  0 - H sin  (p+  #)  - D sin  (p  + t)  cos  at 
S (15) 

+ (Y  + F„)  cos  (P  + t) 

Moment  equilibrium  in  trimmed  flight  is  given  by  EL  = 0,  EM  = 0 which 

s s 

can  be  expanded  using  equations  (13)  and  (l4)  as 

L cos  (p+  | ) - M sin  (p  + f)  + rQ{-  Fa  sin  0 - H sin  (P  + f ) 

- D sin  (p  + * ) cos  a + (Y  + F*)  cos  (p  + * )}  =0 

(16) 


L sin  (P  + Jr)  + M cos  (P  + f)  + rQ  { - F§  sin  0 cos  0 + H cos  (P  + ^) 
+ D cos  (P  + t)  cos  at  + (Y  + F„)  sin  (P  + f)}  = 0 . 

These  equations  may  be  written  in  a somewhat  simpler  form  using 

the  force  equilibrium  conditions  (X_  = 0,  Y = 0)  as 

s s 


X*  - W sin*  0 cos  0 
8 

Y*  * - W sin  0 
s 

Equations  (16)  take  the  somewhat  simpler  form 


(17) 


L cos  (p  + if)  - M sin  (p  + if)  - rQ  W sin  0 = 0 

(18) 

L sin  (p  + if ) + M cos  (p  + if)  - rQ  W sin  8 cos  0 = 0 

Thus  the  trim  conditions  of  the  vehicle  are  given  by  solution  of  equations 
(8)  taken  either  with  equations  (l6)  or  (18).  To  solve  the  trim  problem 
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the  forces  and  moments  appearing  in  these  equations  must  now  be  related 
to  the  flight  velocity  and  control  positions.  The  initial  flight  condition 
is  specified  in  terms  of  horizontal  and  vertical  flight  velocities  with 
respect  to  the  earth,  and  the  gondola  yaw  angle  with  respect  to  this  flight 
path  and  then  the  five  equilibrium  equations  given  above  sire  solved  to 
determine  the  equilibrium  values  of  the  pitch  attitude,  roll  attitude, 
collective  pitch  and  cyclic  pitch  from  these  five  equations.  The  sixth 
equation,  the  yawing  moment  equations  has  not  been  included  since  it  would 
be  used  to  determine  the  power  required.  It  is  assumed  in  the  trim  cal- 
culation that  the  rotor  RIM  is  known. 

k. ) Velocity  Components 

First,  the  velocity  components  must'  be  expresssed  in  various  axis 
systems.  The  flight  condition  is  specified  in  the  gravity  axis  system  by 
the  horizontal  velocity  UQ  and  the  vertical  velocity  WQ,  thus  the  vector 
of  velocity  components  in  the  gravity  axis  system  is 


V0  = 


Uc 

0 

w 


(19) 


The  transformations  given  by  equations  (3)  are  employed  to  find  the 
velocity  components  in  various  axis  systems . 


V,  = V 0 i Va 
V*  - B V, 

Vw  -AV« 

V *8}  V, 
s 


(20) 
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The  sideslip  angle  is  defined  by  the  fact  that  the  velocity  component  v. 


is  equal  to  zero,  i.e., 


V.  - 0 


The  angle-of-attack  is  defined  by  the  fact  that  the  velocity  exponents  vw 


qnH  w are  equal  to  zero,  i.e., 
w 


Vw=  ° 


where  V is  the  resultant  flight  velocity.  From  equation  (20) 


V.  - B * 0 « VQ 


Expanding  equation  (23)  the  following  results  are  obtained  for  the  velocity 
components  in  the  rotor  axis  system,  expressed  in  terns  of  the  horizontal 


and  vertical  velocity, UQ  and  WQ. 


uR  = u0  cos  (P  + t)  cos  e + W0  {-  COS  (P  + ♦)  sin  0 cos  0 
+ sin  (P  + * ) sin  0} 


“o 


- U sin  (p  + t)  cos  9 + W {sin  (p  + * ) sin  9 cos  0 


+ cos  (P  + ♦ ) sin  0} 


w„  = sin  0 + WQ  cos  9 cos  0 


The  sideslip  angle  is  determined  from  the  condition  that  v,  = 0,  so  that 
the  following  relationship  exists 


- Uq  sin  0 + if)  cos  9 + Wq  [sin  (B  + If)  sin  0 cos  0 
+ cos  (B  + f ) sin  0}  = 0 


(25) 


In  the  case  of  level  flight  (WQ  = 0)  this  relationship  simplifies  to 

UQ  sin  (B  + t)  cos  9=0  (26) 

and  therefore  in  level  flight, 

B = - * (27) 

In  general,  equation  (25)  must  be  included  with  the  trim  conditions  to 
determine  the  sideslip  angle.  Once  the  sideslip  angle  is  determined, 
the  angle -of -attack  is  given  by  the  condition 


a = tan-  — (28) 

where  w„  and  u„  are  determined  by  equations  ( 2b ).  Again  in  the  level 

t 

flight  case  (WQ  = 0),  using  equation  (27),  the  following  relationship 
exists 

a = 9 (29)  \ 

i 

The  magnitude  of  the  velocity  in  various  axis  systems  is  given  by  \ 

J 

V = i/ll*  + Woa  =/u,*  + wR8  I 

5. ) Cyclic  Pitch  I 

The  expressions  for  the  rotor  forces  and  moments  are  developed  in  a wind  | 

axis  system  and  consequently  the  cyclic  pitch  included  in  these  equations  is 


referenced  to  the  azimuth  angle  measured  from  the  negative  x„  axis  called 
f ft.  The  physical  system  on  the  vehicle  uses  the  negative  direction  of 
longitudinal  gondola  or  body  axis,  xB,  as  a reference.  This  azimuth  angle 
is  denoted  tB.  Therefore 

te=t*+P  (30) 


The  rotor  or  wind  referenced  cyclic  is  therefore 

AS  = - Al  w cos  t B - B1H  sin  f n (31) 

The  physical  or  shaft  referenced  cyclic  is  therefore 

A9  = - Axb  cos  tB  - Blt  sin  if  B (32) 

Substituting  equation  (30)  and  (31) 

A0  = - [A1M  cos  0 - B1H  sin  0 ] cos  g 

(33) 

- [b1m  cos  0 + A1M  sin  0]  sin  \jrB 


Therefore  the  relationships  between  the  actual  cyclic  pitch  controls  of 
the  vehicle  A1B  and  B10  and  the  wind  referenced  controls  appearing  in  the 


rotor  equations  are 


a1b 

= A1m 

cos  0 

sin  0 

= B1H 

cos  0 

+ Am 

sin  0 

In  particular  these  equations  must  be  noted  when  calculating  stability 
detrivatives . 

6. ) Analytical  Models  For  Forces  and  Moments 

The  expressions  for  the  various  forces  and  moments  described  above 
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must  new  be  expressed  in  terms  of  the  various  flight  condition  variables. 

a. )  Gravity  and  Buoyant  Forces 

The  weight,  W,  is  the  weight  of  the  vehicle  including  the  weight  of  the 
helium  in  the  centerbody.  The  buoyant  force,  Fs,  is  equal  to  the  density  of 
the  air  times  the  displaced  volume 

Fa  = PS  * (35) 

b. )  Rotor  Forces 

A detailed  derivation  of  the  rotor  forces  is  given  in 
Reference  1.  They  are  based  on  the  assumption  that  the  rotor  blades 
do  not  flap  and  can  be  assumed  to  be  infinitely  rigid.  The  rigidity  of 
the  rotor  is  accounted  for  in  the  aerodynamic  model  by  assuming  that 
harmonic  inflow  components  X L and  X H are  developed  proportional  to  the 
aerodynamic  hub  moments  developed  by  the  rotor  such  that 


X m - j ~~~ 
a a 


2Cl  * 

X , = j 

*■  so 


Further  an  harmonic  component  of  the  rotor  inflow  arises  from  the 

"blow  back"  of  the  wake.  The  magnitude  of  this  effect  is  taken  to  be  twice 
that  given  in  Reference  2. 

Xx=-2^tan|  (37) 


where 


-1 

X = tan  x (-4  ) 


Twice  the  value  given  in  Reference  2 was  chosen  for  reasons  discussed  in 
Reference  1. 
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The  total  rotor  inflow  is  given  by 


X =\  - (\x  + X„)  cos  U - Xt  sin  t r 

s 


(39) 


where 


s QR 

and  the  rotor  induced  velocity  is  given  by  momentum  theory  (Ref.  3) 


/,w?  - v 

a V ^ QR  ' + QR 


(40) 


2f 


u. 


The  rotor  advance  ratio  is  u = — . The  various  f factors  account  for  the 

s QR 

fact  that  the  rotor  root  is  at  the  radius  of  the  centerbody,  RB  . 

Define 


B. 

* - T > fn  * 1 - x 


n 


The  expressions  for  the  rotor  forces  in  coefficient  form  are: 


s.!»r 

f fl  J 

\ - *s 

B1H\ 

— - I + 

u f 

a o 3 L 

a f»J  *\ 

1 2 

2 K L ^s  r’ 

r*-1- f>  , .1 

r 3 

1 

a a 2 

L—  w. 

+ -T r 

LXs  * +2 

*3  X«-  flJ 

Aih  f 1 5 

-S-  [*«  fi  +-*aj  ' Xs  f*  + 2i^  ** 


Vi  5 

+ 


rXxfa  ^"1 

Al* 

r . 1 

L 3 1 2 J 

+ ~r 

L*.*  +!Xuf!j 

i«  r 

1 , 

X1  1 

|_XM  fI  + ” 1 

*J  ' Xs 

_x«  fi  +Tf»J 

(hi) 
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c.)  Rotor  Moments 

The  hub  moments  produced  by  the  rotor  were  also  developed  in  Reference 


l 


I 


i 

i 

j 

t 

l 


1,  and  are  given  by 

aih 
2C«  _ ~5~ 
a a 


1 + * f 


k 


1 Sv  _ 

g TTf* 

l + 4 f , 


77*  |--r:tf‘  + il‘3S  f«]  +“a[ff3  + ffJ} 

1 ! T 1 - 5 F f‘ 

+ tttz 


(42) 


d.)  Centerbody  Forces 

The  drag  and  Magnus  forces  acting  on  the  centerbody  are  given  by  the 
expressions 


D = h P (U2  + Wa)  S C, 


(43) 


and 


F„  = | p (Ua  + Wa)  S CL1 


(44) 


The  drag  and  Magnus  force  coefficients  are  assumed  to  be  independent  of 
centerbody  advance  ratio  based  on  limited  data  presented  in  Reference  4. 

This  completes  the  development  required  for  prediction  of  the  trim 
conditions  with  the  exception  of  the  sling  load  and  umbilical  cable  effects. 
These  effects  are  considered  in  a later  section. 

To  summarize  the  trim  calculation,  the  physical  parameters  of  the 


i 


I 

( 

j 


I 
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vehicle  are  specified  and  the  flight  condition  chosen.  The  flight  con- 
dition implies  that  the  horizontal  and  vertical  velocities,  UQ  and  Wq, 
are  chosen  as  well  as  the  body  or  gondola  yaw  angle,  + • The  gondola  yaw 
angle  physically  reflects  the  orientation  of  the  gondola  with  respect  to 
the  trimmed  flight  direction.  Then  the  following  equations  must  be  solved 
simultaneously  in  the  general  case. 

Balance  of  Forces,  equations  (8) 

Balance  of  Moments,  equations  (16)  or.  (18) 

Velocity  Component  Relationships,  equations  (24) 
with  the  conditions 

VR  - 0 

_1  WR 
a = tan  — 

Vehicle  control,  wind  referenced  control  relationship, 
equations  (34) 

Rotor  and  Centerbody  Forces  and  Moments  in  terms  of 
Flight  Condition,  equations  (4l),  (42),  (43),  and  (44) 

The  solution  of  these  equations  yields  trim  values  of  the  rotor  cyclic 
and  collective  pitch,  the  vehicle  pitch  and  roll  angles,  and  the  angle-of- 
attack  and  sideslip  angle. 

The  velocity  relationships  are  considerably  simplified  in  the  case 
of  level  flight  reducing  to  the  conditions  0'=-t,  and  or  = 0. 

It  should  be  noted  that  this  equilibrium  solution  does  not  include 
the  effects  of  the  sling  load  and  umbilical  cable.  Appendix  D discusses 
incorporation  of  these  effects  into  equations  (8)  and  (16). 
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Level  flight  was  characteristic  of  all  the  flight  conditions 
examined  in  the  experimental  program.  Further,  in  the  flight  program, 
the  gondola  orientation  was  always  selected  such  that  f = 0 and  therefore 
from  (27),  9 ■ 0. 

Consequently  the  simplified  trim  equations  for  initially  level  flight 
are,  placing  condition  (27)  on  equations  (8)  and  (16)  level  flight  equations 
are  as  follows: 

(-W  + Fb)  sin  0 cos  0 - H - D cos  8=0 
(W  - F,)  sin  0 + (Y  + F„)  = 0 

(W  - FB)  cos  0 'cos  0 - T - D sin  0 = 0 (45) 

L + rQ{-  F,  sin  0 + (Y  + F„)}  = 0 
M + rQ{-  FB  sin  0 cos  0 + H + D cos  0}-=  0 

The  last  two  equations  may  also  be  written  alternatively  using  equations  (18) 
as 

L - r W sin  0 = 0 

° * (U6) 

M - r W sin  0 cos  0=0 
o 

These  equations  are  then  solved  using  the  rotor  and  centerbody  force  and 
moment  relationship  given  by  equations  (4l)  and  (42).  Note  also  that  since 
the  initial  sideslip  is  zero 

Ais  * Aiw  (47) 

= Bit* 

The  equations  given  by  (45)  were  programmed  on  a digital  computer  and 
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solved  for  the  trim  condition.  The  method  of  approach  was  to  initially 
solve  the  three  force  equations  neglecting  the  rotor  in-plane  and  side 
forces  to  determine  initial  values  of  the  roll  and  pitch  attitude  and 
then  the  moment  equations  can  be  solved  for  the  cyclic  pitch  required. 

Once  an  initial  trim  solution  is  obtained  the  rotor  in-plane  and  side 
forces  are  calculated,  added  to  the  equations  and  a new  trim  condition 
calculated.  Typically  this  procedure  converges  rapidly  as  the  rotor 
in-plane  and  side  forces  have  only  a small  influence  on  the  trim  attitudes. 
II.  DYNAMIC  STABILITY  ANALYSIS 

This  section  develops  equations  of  motion  for  the  AEROCRAME  incor- 
porating a sling  load.  The  attachment  point  of  the  load  on  the  vehicle 
is  taken  to  be  an  arbitrary  point  located  on  the  axis  of  rotation.  The 
load  is  assumed  to  be  a point  mass  and  has  two-degrees -of-freedcm  with 
respect  to  the  body.  Aerodynamic  forces  on  the  sling  load  are  neglected. 

The  forces  and  moments  acting  on  the  vehicle  were  resolved  into  a 

shaft  axis  system(X  *)  in  section  I.  The  shaft  system  orientation  is 
s 

obtained  by  rotating  the  gravity  axis  system  (XQ)  through  the  vehicle  roll 
and  pitch  angles.  We  now  introduce  two  additional  axis  systems  as  necessary 
to  proceed  with  the  development  of  the  vehicle-sling  load  equations: 

a.  Attachment  Point  Axes  X*  (x*,  y*,  zA) 

These  axes  are  always  parallel  to  the  shaft  axes  X . The  origin  of 

s 

this  axis  system  is  located  at  the  attachment  point. 

b.  Load  Axes  Xt  (xL,  yL,  zL) 

This  axis  system  is  a body  axis  system  attached  to  the  load  with  its 
origin  at  the  load  center  of  gravity.  The  zL  axis  points  downward  and  is 
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parallel  to  the  support  cable.  These  axes  are  shewn  schematically  In  Figure 

1(c)  and  Figure  1(d)  shows  the  freebody  diagram  employed  In  the  analysis. 

X and  L are  the  force  and  moment  vectors  acting  on  the  vehicle  form- 
s s 

ulated  in  Part  I given  by  equations  (7)  and  (10).  Rt  is  the  force  vector 
which  represents  the  force  applied  to  the  vehicle  produced  by  the  load  and 
acts  at  the  attachment  point. 

The  equations  of  motion  for  the  vehicle-load  system  can  be  written  in 
vector  form  using  Newton’s  Laws  as 

Hco|r  +^s  x HC0|f=  PA  x + Lg 

»'  aca|,  = Rt  + X.  (48) 

mi  . ®c ajL  = - RL  + Rsl 
Hco|l  + H L x Hcq|i.=  (“  x (“  Rl) 

The  first  two  equations  are  the  equations  of  motion  of  the  vehicle  where 

X and  L were  calculated  in  a previous  section.  The  third  and  fourth 
s s 

equations  are  the  load 'equations  of  motion.  It  has  been  assumed  that  there 
are  no  external  forces  acting  on  the  load  with  the  exception  of  the  gravity 
forces.  This  assumption  is  quite  reasonable  owing  to  the  high  density  of 
the  load  employed  in  the  experiments.  The  sign  of  the  moment  term  in  the 
last  equation  arises  from  the  fact  that  moments  are  taken  about  the  load 
center  of  gravity.  Hc«|r  is  the  moment  of  momentum  of  the  vehicle  with 
respect  to  the  center  of  gravity  of  the  vehicle  and  is  the  angular  rate 
of  the  Xs  axis  system.  Hc 9|L  is  the  moment  of  momentum  of  the  load  with 
respect  to  the  load  center  of  gravity  and  L is  the  angular  velocity  of  the 
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load  axis  system  with  respect  to  space.  The  support  cable  has  been  assumed 
to  be  a massless  rod.  It  is  further  assumed  that  the  load  is  a point  mass 
that  is,  that  its  moments  of  inertia  are  zero  with  respect  to  its  center 
of  gravity  and  consequently  HCbjl  = 0.  Therefore  the  fourth  of  equations 
(48)  becomes 

PL  x RL  = 0 

This  implies  that  the  reaction  force  RL  always  lies  along  the  support 
cable  direction.  The  reaction  force  can  therefore  be  eliminated  from  the 
third  of  equations  (48)  by  taking  the  cross  product  with  PL  so  that  the 
third  equation  becomes 

mL  (PL  x aCajL  ) = Pt  x FaL 

RL  is  then  eliminated  from  the  first  and  second  equations  using  the  third 
equation 

= Fau  - mL  aCQjL 

so  that  equations  (48)  can  be  written  as 

Hcajr  + x ^cq|f  * P*  x ” ®i.  aco|L  ) + 

aca|r  » (faL  - mL  acajt  ) + Xg  (49) 


® L X ®C  aju  ) = X ^®L 

the  last  equation  can  also  be  written  as 

PL  x (mL  aCajL  - FaJ  = 0 (49a) 

Now  the  individual  terms  in  equations  (49)  and  (49a)  are  developed.  I,  j , 
k denote  unit  vectors  with  subscripts  indicating  their  axis  orientation. 
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The  various  accelerations  involved  in  the  equations  of  motion  may  be 


u= 


cos  0 sin  9 sin  0 - sin  0 cos  0 


cos  0 


sin  0 


sin  9 - sin  0 cos  0 cos  6 cos  0 


cos  0 1 sin  0 L sin  0 L - sin  0 L cos  0 , 


cos  0L 


sin  0L 


sin  0 L - sin  0 L cos  0 L cos  0 L cos  0 L J Ik 


The  development  is  specialized  for  hovering  flight  by  assuming  that 
the  initial  vehicle  and  load  attitude  angles  are  zero  such  that  0 , 0 , 0 L 
and  0 L represent  perturbations  from  trim.  It  is  further  assumed  that 
these  perturbations  are  small  so  that  equations  (52)  can  be  expressed  as 


is  - i.  - 9 K 


3S  - 30  + <t>  ka 


kg  = 9 i a - 0da  + ka 
- 13  - «t  t, 
h - ia  * *L  % 


kL  = 0L  is  - 3S  + *a 


For  the  load  equation,  it  is  desirable  to  express  the  gravity  force  in 
terms  of  components  in  the  load  axis  directions . This  is  accomplished 


J 


by  inverting  transformation  (52) 


cos  8 


sin  6 


sin  9 sin  0 cos  0 - sin  0 cos  9 

-sin  9 cos  0 sin  0 cos  9 cos  0 


cos  9l 


k 

sin  9l  sin  0L  cos  0^ 


sin  0t 

sin  0L  cos  0L 


i 

-sin  0L  cos  0l  sin  0L  cos  0L  cos  0^ 


For  small  angles 

- K - 8 £„ 

= -<»ss 

£«  - - e h * 4 *.  * K (5, 

I =.  IL  + 0L  kL 

S 9 

3S  - 3t  - \ 

£s  - - *<■  it  + *.  3t  + £i 

Substituting  from  tbe  second  set  of  relationships  in  (5*0  into  the  first 


£a  = ±L  (-9  - 9L  ) + (0  + 0i  ) 5t  + £l 


where  only  first  order  terms  have  been  retained  as  consistent  with  the  small 
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angle  approximation. 

The  load  angular  velocity  expressed  in  load  coordinates  becomes 

fit  = (Ps  + Pt)  it  + (q.s  + <li)  *t 

where  only  first  order  terms  have  been  retained. 

The  load  acceleration  is  given  by 

aca|L  = ®a  + (4S  + 4l)  Zl  it  “ (Ps  + Pt)  Zt  Jt  (56) 

where  again  only  first  order  terms  have  been  retained.  The  acceleration 
at  the  attachment  point  aA  is  given  by 


at  ~ aca  r + f4flZA  ig} 


(57) 


- fisz*  3.) 


Converting  equation  (57)  to  the  load  axes  reference  using  (54),  retaining 
only  first  order  terms  and  incorporating  into  (56), 


a c g| t = acs|r  + (4S  (ZA  + ZL)  + qL  ZL]  IL 
” (Pg  (Z*  + Zj  + PL  Zj 

Equations  (58)  and  (55)  are  now  employed  to  express  the  reaction  force 
Rt  in  load  coordinates 

™ ^Gt  “ “-i  ac ajt 


(58) 


Rl  = mL  g {kL  - (9  + eL ) IL  + (0  + 0L)  3J 


- {aC8|f+  [4g  (ZA  + Z J + 4l  zj  IL  (59) 

- C*8  (ZA  + zj  + ZJ  i J 


Expressing  the  components  of  this  vector  as 

XL  = - mL  g (8  + 9t  ) - mt  [q  (ZA  + Zj  + 4t  Z J - mL  u 

s s 

YL  = *L  g (0  + 0,.  ) + aL  [pg  (ZA  + Zj  + pL  Z J - mL  vg  (59a) 


Zu  = mL  g 
so  that 

(^GL  " ac  s|l  ) = i Y^  j ^ + Z k|_ 

The  load  equation  of  motion  is  obtained  from  (49a)  as 
P]_  x (t#l  aCQjL  “ Fet  ) = 0 


or 


3l  Xl  - IL  ZL  Yt  =0 

For  the  second  of  equations  (49)  the  reaction  Rt  is  expressed  in  terms  of 
shaft  coordinates.  This  is  obtained  from  the  unit  vector  relationships 
given  in  (53)  taken  with  equations  (59a).  Since  Xu  and  YL  are  perturbation 
quantities  the  result  is 

(Fat  - mt  aCQjL  ) ■ (XL  + 9t  Zj  Is 


(60) 


+ (Yt  - 0l  Zu)  3a  + ZL  kg 
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since  ZL  = mLg,  equation  (60)  can  be  written  as 
Fgl  - ®u  aca|L  = (Xt  + aLg  9l  ) ig 


(60a) 


+ (Yl  - aLg  0U)  3S  + mLg  kg 


This  is  the  form  of  the  reaction  force  required  for  the  second  of 
equations  (49).  For  the  first  of  equations  (49),  the  cross  product  PA  x 
(Fgl  - mL  ac„jL  ) must  be  calculated  where  PA  = zA  kg.  Taking  the  cross 
product  using  equations  (60a) 

Pa  x (F9L  - mL  aCa|L  ) = 3S  (ZA  (Xw  + mLg  0L  )} 

- ig  (ZA  (Yl  - mLg  0L  )} 

The  equations  of  motion  can  be  written  as 


Hca|r  + ns  x Hca|,  = 3g  £ZA  (Xu  + mLg  0L  )} 


- ia  [Z*  (Yu  - mug  0l  )}  + Lg 


aCa(r  * (xl  + °HS  9u  ) 


+ (Yl  - toLg  ) 3g  + mug  kg  + Xg 


A L Zl  Xu  - it  Zu  Yu  = 0 

where  XL  and  Yt  are  given  by  equations  (59a). 

Equations  (6l)  are  the  equations  of  motion  for  the  AEROCRANE  with 
a point  mass  sling  load.  These  equations  can  be  simplified  by  noting 
that  the  last  of  these  equations  reduces  to 
X,  = 0 

( 

Yu  = 0 
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so  that  the  first  two  equations  become 

Hca|f  + x Hc 0| r = Ls  + jg  (ZA  mLg  0t  ) 

+ Ig  (ZA  mLg  0L ) 

m.'  acajf  = mLg  0L  ig  - mLg  0l  3g  + *\JB  + Xg 

Expanding  equations  (62)  and  (63) 
l'p  + \ Oq  = Lg  + ZA  mLg  0L 

I#4  - Iz  Qp  - Mg  + ZA  mLg  0L 

»'  % - Xg  + »Lg  0, 

, ( 

“ *B  = YS  " “Lg  0 L 

- mL  ug  - mLg  (9  + 0L  ) - mL  (4g  (ZA  + ZL)  + 4L  ZL)  = 0 

- mL  vg  + mLg  (0  + 0L  ) + mL  (pg  (ZA  + Zt)  + pL  ZL)  = 0 

Equstions  (64)  are  equations  of  motion  for  the  AEROCRANE  including  the 

# 

effects  of  a sling  load.  They  have  been  specialized  for  the  hovering 
case  by  the  assumption  of  zero  initial  attitudes.  For  the  general  case 
of  forward  flight  the  assumption  of  zero  initial  trim  angles  must  be 
removed. 

as  a result  of  the  small  angle  assumption 


i « 


(65) 


Substituting  into  equations  (64) 

i'  0 + I Q 0 = L + ZA  mtg 
z s 

i'  0 - I n 0 = M + ZA  mLg  0L 

Z S 

a'  &s  = Xg  + mLg  0L 
«'  = Ys  - mLg  <K 

- mL  us  - mug  (9  + 0L  ) - »L  (0  (ZA  + Zu)  + 0L  ZL)  * 0 

- mL  V + mLg  (0  + 0L  ) + aL  (0  (ZA  + Zt)  + 9^  ZL)  = 0 

The  expansion  of  X„ , Y , L , and.  M in  a Taylor  series  about  the 

s s s s 

equilibrium  flight  condition  has  been  treated  in  detail  in  Reference  1. 
There  are  additional  aerodynamic  terms  which  were  developed  in  Reference 
1 which  depend  upon  the  acceleration  of  the  vehicle  and  arise  from  the 
fact  that  the  center  of  gravity  of  the  vehicle  is  not  located  at  the 
center  of  buoyancy. 

These  terms  are  ^ith  signs  for  the  right  hand  side  of  the  equations 

Rolling  Moment  equation 

- r m v\ 
o a s 

Pitching  Moment  equation 


r„  K 

o a s 

Horizontal  Force  equation 

r « 0 

o a 

Lateral  Force  equation 


r m 0 
o a 


(66) 
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Further  the  effects  of  the  umbilical  cable  must  be  added.  These  are 
developed  in  detail  in  Appendix  B and  are  given  by  the  sum  of  equations 
(B-l8)  and  (B-3*0  as 

Rolling  Moment  equation 

<$ 

“c  z*  <-*s  + z*  •)  - V*  + » 

Pitching  Moment  equation 

d0 

"ce  *.  (4.  + Z.  «)  ' V*  (1  + df>  8 

Horizontal  Force  equation  (67) 

”c.  <V+Z> s>  - wc  (1  + dr) 6 
Lateral  Force  equation 

* 

mce  (*s  “ Z*  ^ + Wc  (l  + d"?  0 

Combining  equations  (65),  (66)  and  (67)  gives  the  hovering  equations  of 
motion. 

These  equations  may  be  written  compactly  in  matrix  notation  as 

M {q}  + C {fl}  + K [ql  « F {5}  (68) 

The  mass,  damping,  spring,  and  control  matrices,  M,  C,  K,  and  F are 
given  on  the  following  pages.  The  motion  variables  [q]  are 

10 

3 

0t 

eL 

x 

y 


The  control  matrix  {8}  is 


These  equations  of  motion  were  solved  to  determine  the  characteristic 
dynamics  of  the  vehicle. 

The  moment  equations  are  normalized  by  the  inertia  I/ , the  force 
equations  by  the  mass  m , and  the  sling  load  equations  are  divided  by 
the  sling  load  mass  times  the  sling  load  length.  In  order  to  make  the 
notation  more  compact  the  following  symbols  were  introduced. 

PB  r 
a _ o 


I 

= — Q 


u)  a _ _g_ 
sl  Zu 

In  addition,  the  symnetry  properties  of  the  aerodynamic  derivatives  were 

employed  to  eliminate  the  lateral  derivatives,  i.e., 

L * M Lb  Ma 

P <1  13  is 


The  subscript  notation  for  the  aerodynamic  derivatives  implies-  that 
the  moment  derivatives  have  been  divided  by  I7  and  the  force  derivatives 


by  m . 
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The  object  of  the  experimental  program  was  to  determine  the  vehicle 
transfer  functions,  which  can  be  obtained  from  equations  (68)  in  the  following 
fashion.  First  the  Laplace  transform  of  equations  (68)  is  taken  giving 

[Ms3  + Cs  + K]  {q3  = F {A}  (69) 

where  [q]  is  the  Laplace  transform  of  the  motion  variables  {q}  and  {a}  is 
the  Laplace  transform  of  the  control  variables  {6} . Equations  (69)  become 

{q3  * [Ms8  + Cs  + K]*1  F {A}  (70) 

There  are  twelve  vehicle  transfer  functions,  given  by  the  elements  of  the 
matrix  [Ms3  + Cs  + K]-^  F.  The  twelve  elements  of  the  matrix  are  shown 
symbolically  in  the  next  page.  These  elements  characterize  the  response 
of  the  vehicle  and  sling  load  to  the  control  inputs.  Since  there  are  3ix 
motion  variables  and  two  control  variables,  twelve  transfer  functions  are 
required  to  characterize  the  dynamic  motions  of  the  vehicle. 

The  transient  response  experiments  were  designed  to  verify  the  analytically 
calculated  transfer  functions  which  describe  mathematically  the  dynamic  response 
of  the  model.  Owing  to  the  polar  symmetry  of  the  model  and  its  aerodynamic 
characteristics  in  hovering  flight,  six  of  these  transfer  functions  are  simply 
related  to  the  other  six  by  the  following  relationships 


££s) 

Ai.(s) 

= 

Bia(s ) 

1(5,  L- 

=.^8iT 

Blt(s) 

Ais(s) 

(s) 

1 (s) 

Aib(s) 

B1#(s) 

1 (s) 

0L  (s) 

A1#(s) 

Bi8(s  ) 
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m 


•‘isl-'  — 

Aii(*)  Bn(a) 


Owing  to  these  relationships  only  a Blt  input  was  applied  in  hovering 
as  the  response  of  the  vehicle  to  An  inputs  can  be  determined  front  the 
response  to  61S  inputs. 

The  experimental  verification  of  these  transfer  functions  is 
discussed  In  a later  section  of  the  report. 
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II.  MODEL 

The  model  employed  for  the  experiments  described  in  this  report  is  a 
! modification  of  the  0.107  Froude  scale  model  of  a proposed  50  ton  payload 

AEROCRANE  vehicle  as  described  in  Reference  1.  The  modifications  consisted 
j principally  of  an  increase  in  the  spherical  centerbody  diameter,  provision 

i of  a dumpable  water  ballast  package  and  incorporation  of  a stability  aug- 

j mentation  system.  According  to  the  conclusions  and  recommendations  made 

j ‘ 

in  Reference  1,  these  modifications  were  provided  to  increase  flight  safety 
by  allowing  positive  buoyancy  to  be  achieved  at  all  times  and  to  ease  the 
remote  pilot's  task  in  controlling  the  lightly  damped  model  motions  with 
inadequate  motion  cues. 

A photograph  of  the  modified  model  in  hovering  flight  and  showing  the 
dumpable  water  ballast  sling  load  is  shown  in  Figure  2 , and  a 2-view  drawing 
. is  presented  in  Figure  3.  Table  I presents  a summary  of  the  model  physical 

characteristics.  « 

IT.  CENTERBODY  MODIFICATIONS 

To  achieve  the  additional  buoyancy  required  for  flight  safety  consider- 
ations the  gas -containing  spherical  centerbody  was  increased  in  diameter  from 
16  feet  to  18  feet,  providing  an  additional  60  lb  of  buoyant  lift.  The 
increased  gas  envelope  size  necessitated  a revised  internal  structure  con- 
| sisting  of  longer  frame  members,  addition  of  internal  stress-relief  cables 

and  gas  envelope  stress-relief  patches  at  the  radial  pass-through  fittings 
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to  accommodate  the  increased  buoyant  gas  loads.  These  structural  modifications 
are  shewn  schematically  in  Figure  4. 

HI.  WATER  BALLAST  SYSTEM 

The  jettisonable  water  ballast  used  to  simulate  a payload  was  carried 
in  a special  aluminum  container  with  a trap-door  type  bottom  and  suspended 
from  the  model  gondola  as  a sling  load  as  shown  in  Figure  2.  The  water  was 
contained  in  a plastic  bag  liner  fitting  within  the  aluminum  container.  The 
trapdoor  latch  was  secured  with  light  polyester  fishing  line  which  was  in 
turn  wrapped  with  a coil  of  high  electrical  resistance  wire.  Upon  command 
from  the  truck-based  operator  an  electric  current  would  heat  the  high 
resistance  wire,  melt  the  polyester  fishing  line  and  allow  the  trapdoor  to 
open,  thereby  jettisoning  the  water  ballast.  Electric  power  for  the  jettison 
was  provided  by  a dedicated  24  volt  battery,  thereby  assuring  operation  even 
with  complete  power  loss  from  the  truck  generator  systems. 

IV.  STABILITY  AUGMENTATION  SYSTEM 

t 

As  a result  of  the  flight  test  experience  and  analysis  efforts  reported 
in  Reference  l.} it  was  determined  that  a stability  augmentation  system  (SAS)  to 
stabilize  the  lightly-damped  retrograde  precessional  mode  would  greatly  ease 
the  remote  pilot's  burden  in  controlling  the  model.  Although  it  was  demon- 
strated in  analog  simulator  flights  that  the  pilot,  with  adequate  motion  cues, 
could  stabilize  this  mode,  in  the  model  flight  operations  the  required  cues 
were  not  available  to  the  remote  pilot,  and  controlled  flight  was  nearly 
impossible.  Further,  it  was  shown  in  Reference  lthat  a phased  attitude 
feedback  given  by  the  expressions 
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A1B  = Ka  [0  sin  y - 6 cos  y] 


(69) 


Bis  * ka  f#  cos  Y + ® sin  y]  (70) 

would  effectively  stabilize  the  lightly-damped  precessional  mode.  Accordingly, 
for  the  experiments  reported  herein,  a stability  augmentation  system  was 
implemented  to  accomplish  this  task. 

In  view  of  the  exploratory  nature  of  the  research,  the  time  constraints 
associated  with  the  experimental  operations  and  the  test  objective  of  forward 
flight  experiments,  a more  general  stability  augmentation  system  was  designed 
and  installed  in  the  model  controller.  This  stability  augmentation  system, 
schematic  representations  of  which  are  shown  in  Figures  5 and  6,  allowed 
for  selection  of  any  desired  feedback  phasing  angle,  y,  in  addition 
to  y = 45°  as  indicated  by  equations  (69)  and  (70).  Additionally, 
although  the  hovering  analysis  of  Reference  1 indicated  no  specific  need 
for  angular  rate  feedbacks,  these  were  allowed  for  in  the  implementation 
as  diagrammed  in  Figure  6. 

9 

Model  attitude  and  rate  information  was  supplied  to  the  SAS  (stability 
augmentation  system)  by  the  three-axis  integrating  rate  gyro  package  des- 
cribed in  Reference  1.  This  instrumentation  package,  originally  intended 
for  flight  dynamic  data  acquisition,  was  more  than  adequate  for  SAS  inputs 
and  performed  faultlessly.  The  analog  integrators  required  for  determining 
model  attitude  from  the  rate  gyro  information  were  revised  to  virtually 
eliminate  environmentally-produced  drift  errors. 
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V.  AZIMUTH-HOLD  LOOP 


The  azimuth- hold  loop,  driving  the  retrograde  motor  that  allowed  the 
gondola  to  be  positioned  azimuthally,  was  revised  to  provide  greater  torque 
capability  and  thereby  eliminate  the  difficulties  experienced  in  the  flight 
tests  report  in  Reference  1.  The  revised  system  possessed  approximately 
4 times  the  torque  capability  of  the  original  system.  Extensive  laboratory 
testing  of  this  vital  loop  closure  was  performed  on  a specifically-designed 
flight  simulation  set-up  to  insure  satisfactory  performance,  and,  although 
in-flight  performance  was  adequate  throughout  the  flight  envelope,  certain 
dynamic  problems  were  encountered  in  flight,  and  loop  compensation  adjustments 
were  required. 

VI.  RAPID  DECELERATION  SYSTEM 

A necessary  function  in  the  flight  safety  systems  that  includes  rapid 
achievement  of  a fully-buoyant  state  through  ballast  dump  is  the  rapid 
deceleration  of  the  model  rotational  motion.  The  model  rpm  control  system 
operates  through  varying  the  speed  of  the  four  fixed-pitch  propellers  mounted 
on  the  model  wings  by  adjusting  the  output  of  the  main  400  Hz  model-power 
alternator.  In  order  to  decelerate  the  model  rotational  motion  rapidly  it 
is  necessary  to  actually  reverse  the  propulsive  motor  voltage  polarity  and 
hence  reverse  the  direction  of  rotation  of  the  propellers.  To  accomplish 
this,  reversing-current  relays  were  installed  on  the  model  gondola  electrically 
downstream  of  the  rectifier  package.  A high  sensitivity  alternator  field 
control  potentiometer  was  incorporated  to  allow  the  operator  to  rapidly 
reduce  the  motor  power  to  near  zero,  actuate  the  polarity  reversing  relays 


t . 

f ' 
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and  subsequently  reactivate  the  alternator  field  to  provide  reversed  thrust. 
The  system  in  operation  was  capable  of  arresting  the  model  rotation  in 
approximately  one  revolution  without  exceeding  model  motor  rated  currents. 
Thus, in  the  event  of  an  in-flight  problem  that  threatened  the  flight 
safety  of  the  model,  the  model  rotational  motion  could  be  arrested  and  a 
fully-buoyant  state  achieved  by  ballast  dump  in  approximately  3 seconds. 


45 


EXPERIMENTAL  FLIGHT  TEST  PROGRAM 


The  experimental  flight  test  program  was  conducted  in  Hangar  No.  1 at 
the  Naval  Air  Engineering  Center,  Lakehurst,  N.  J.  in  the  time  period  from 
March  21,  1977  to  May  2,  1977.  A total  of  56  data  runs  were  performed  and 
approximately  20  hours  of  flight  time  were  accumulated.  The  principal 
test  objectives  were  to  quantify  the  model  trim  conditions  and  transfer 
functions  as  they  varied  with  flight  condition,  model  configuration  and 
stability  augmentation  (SAS)  for  correlation  with  the  theoretically  pre- 
dicted characteristics. 

The  flight  testing  efforts  and  procedures  were  divided  into  two 
principal  types  of  tests,  hovering  and  forward  flight;  each  of  these  types 
will  be  discussed  separately. 

I.  HOVERING  FLIGHT 

A photograph  of  the  model  in  a typical  hovering  flight  is  presented  in 
Figure  7 showing  the  model,  water  ballast  sling  load  and  umbilical.  At  the 

f 

bottom  of  the  photograph  can  be  seen  the  parked  truck  which  carries  the  crew, 
model  control  and  power  systems.  A photograph  of  the  flight  crew  arrangement 
is  presented  in  Figure  8 showing  the  pilot,  flight  engineer  and  test  director 
at  the  model  system  control  and  data  consoles.  Although  the  photograph  .of 
Figure  7 shows  the  bottom  of  the  umbilical  suspended  from  a tower  mounted 
on  the  test  truck,  most  of  the  hover  runs  were  accomplished  without  the 
tower , and  the  umbilical  ran  downward  from  the  model  directly  to  the 
ground.  This  arrangement  allowed  hovering  tests  to  be  performed  without 
any  horizontal  force  or  moment  contribution  from  the  umbilical  cable 


Other  hover  tests  with  the  umbilical  tower,  as  pictured  in  Figure  7,  were 
performed  with  horizontal  force  and  moment  initial  conditions  imposed  by 
the  umbilical  catenary  shape.  a summary  of  the  hovering  test  conditions 
investigated  is  presented  in  Table  II. 

The  hovering  tests  were  performed  by  establishing  a steady  hover 
condition  with  the  bottom  of  the  model  approximately  90  ft.  above  the 
ground  and  with  the  umbilical  hanging  directly  downward.  This  zero- 
initial-condition  hover  was  established  with  the  SAS  in  operation  at 
y = 45°  and  with  the  gain  KA  = 0.3  °/o*  This  SAS  configuration  has  been 
established  on  the  initial  run  as  being  a level  of  stability  augmen- 
tation that  was  completely  comfortable  to  the  pilot  and  for  which  the  un- 
disturbed model  motions  were  undetectable.  At  these  initial  conditions  the 
SAS  gain  was  reduced  to  the  level  desired  for  the  particular  test  sequence, 
and  a pulse  input  in  cyclic  pitch  was  applied  by  the  test  engineer  by 
means  of  a switch  on  the  engineer's  console.  After  establishing  the  initial 
conditions  and  throughout  the  ensuing  transient  response,  the  pilot's  controls 
were  held  fixed.  Figure  9 shows  typical  hovering  transient  responses  for  the 
unstabilized  model  and  with  various  level  of  SAS  gain. . 

II.  FORWARD  FLIGHT  EXPERIMENTS 

Photographs  of  the  model  during  forward  flight  experiments  are  presented 
in  Figures  10  and  11.  In  these  figures  can  be  seen  the  relative  position  of 
the  model,  umbilical  catenary  and  the  truck  carrying  the  crew  and  model  support 
systems.  The  forward  flight  tests  required  additional  instrumentation,  not 
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necessary  for  the  hovering  experiments,  consisting  of  a fifth  wheel  for 
measuring  truck  speed,  and  umbilical  shape  and  position  instrumentation. 

The  fifth  wheel  is  mounted  under  the  rear  of  the  truck  as  seen  in  Figure 
11  and  the  umbilical  shape  and  position  instrumentation  is  mounted  at  the 
top  of  the  umbilical  support  tower  on  the  truck. 

The  test  procedure  for  the  forward  flight  runs  commenced  with  es- 
tablishing a steady  hover  in  the  northwest  corner  of  the  hangar  immediately 
to  the  right  of  the  center  crack  in  the  hangar  doors  visible  in  Figure  10. 
Owing  to  the  presence  of  a fenced  in  storage  area,  also  visible  in  the 
foreground  of  Figure  10,  it  was  necessary  to  fly  along  a diagonal  path 
towards  the  southeast  corner  of  the  hangar.  Runs  were  made  in  one  direction 
only,  principally  to  eliminate  the  necessity  of  reccmpensating  the  rate 
gyro  package  for  earth's  rotational  rate.  Once  a steady  hover  had  been 
established  with  the  model  in  a position  relative  to  the  truck  that  was 
acceptable  with  respect  to  altitude  and  umbilical  shape,  the  pilot  would 
apply  the  cyclic  and  collective  pitch  inputs  required  to  transition  the  model 
to  the  desired  forward  flight  condition.  The  truck  driver  was  required  to 
adjust  the  truck  speed  to  stationkeep  with  the  model.  The  desired  forward 
flight  trim  conditions  were  established  with  a nominal  SAS  gain  of  approxi- 
mately 0.3  °/  and  y = 45°.  For  the  initial  forward  flight  runs  the  trim 
condition  was  maintained  steadily  with  no  additional  control  inputs  in 
order  to  measure  the  model  trim  conditions.  At  a predetermined  location 
along  the  line  of  flight  the  transition  back  to  hovering  flight  was 
instituted  and  the  end  of  run  hover  established.  In  general,  excepting 
for  the  very  lowest  speed  flights,  transitions  to  and  from  the  forward 
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flight  trim  conditions  were  made  with  the  pilot  operating  in  an  open  loop 
fashion.  That  is,  there  was  inadequate  space  within  the  hangar  for  the 
pilot  to  perform  a transition  by  a series  of  small  perturbations  to  the 
controls  and  subsequent  corrections  to  flight  path  errors.  Instead,  estimated 
trim  control  positions  were  predetermined  and  the  pilot  simply  put  the 
control  at  these  values,  at  a prudent  rate,  to  perform  the  transition.  The 
SAS -augmented  model  stability  level  was  such  that  no  difficulty  was  experi- 
enced in  transitioning  in  this  manner.  A typical  transition  time  history 
is  presented  in  Figure  12  shewing  the  pilot 1 s control  Inputs  and  the  ensuing 
model  pitch  and  roll  Angni motions. 

The  data  required  to  determine  the  model  transfer  function  to  control 
inputs  were  obtained,  as  in  hover,  by  applying  an  electrical  pulse  input 
to  the  cyclic  controls  by  means  of  a switch.  In  forward  flight,  however, 
in  order  to  obtain  the  unstabilized  (K*  = 0)  model  transfer  function  and 
still  maintain  the  required  trim  control  settings,  it  was  necessary  to  utilize 
track  and  hold  networks  on  the  model  control  signals.  These  networks  held 
the  trim  controls  required  while  the  SAS  gain  was  switched  to  zero  to  obtain 
the  desired  transient  response  time  histories.  Such  a typical  transient 
response  time  history  at  a forward  flight  trim  condition  is  shown  in  Figure  13. 
Table  III  present  a summary  of  the  forward  flight  trim  conditions. 


where  the  additional  terms  SX  and  SY  represent  the  measured  umbilical  forces 
acting  on  the  model, and  level  flight  at  zero  slideslip  has  been  assumed. 
These  quantities  have  been  evaluated  using  the  experimental 3 y-  measured 
values  for  all  terms  excepting  Y and  H,  the  rotor  in-plane  forces,  which 
were  evaluated  using  the  theoretical  model  at  the  experimental  operating 
conditions.  The  resulting  values  are  presented  in  Figures  l4  and  15  and 
compared  with  the  theoretical  values  for  FM  and  D employing  the  assumed 
force  coefficients. 

As  can  be  seen  in  Figure  l4,  the  original  value  of  = 0.6  seriously 
under  predicts  the  values  of  centerbody  drag  extracted  from  the  experimental 
data  and  a value  of  CD  = 0.80  is  chosen  as  being  a fair  representation  of 
the  drag  coefficient  within  the  experimental  scatter.  The  sources  of  the 
experimental  scatter  in  the  drag  data,  include,  in  probable  order  of 
importance,  determination  of  model  velocity  from  truck  speed  measurements, 
wind  currents  within  the  hangar  and  accuracy  of  measurement  of  umbilical 
shape  and  position. 

# 

The  last  source  of  error  is  particularly  important  in  the  determination 

of  F„  since  the  umbilical  lateral  force  term,  SY,  was  not  measured  directly 

but  was  inferred  from  the  beginning-of-run  hover  umbilical  position.  This 

uncertainty  is  reflected  in  the  relatively-larger  scatter  shown  in  the 

experimental  F*  data  presented  in  Figure  15.  The  theoretical  Magnus  force, 

. ( assuming  a value  of  C-  * 0.30,  is  shewn  in  Figure  15  in  comparison  with 

[ ' T4 

the  experimental  data.  Due  to  the  mentioned  experimental  uncertainties  and 

the  resulting  data  scatter  it  was  concluded  that  it  would  be  unjustified  to 

attempt  to  refine  the  Magnus  lift  coefficient  value  any  further. 

I ■ 
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2.  Trim  Control  Requirements 

The  data  presented  In  Figure  16  show  comparisons  of  theoretical  and 
experimental  values  for  the  longitudinal  and  lateral  cyclic  pitch  required 
for  trim.  While  the  force  equilibrium  equations  (8)  used  to  express  the 
centerbody  aerodynamic  forces  are  influenced  by  the  assumed  rotor  aero- 
dynamics through  the  H and  Y terms,  the  cyclic  pitch  requirements  for  trim 
presented  in  Figure  16  are  much  more  indicative  of  the  accuracy  of  the 
rotor  aerodynamics  representation.  In  particular  the  good  low  speed  agree- 
ment between  theory  and  experiment  shown  in  the  A1S  data  indicates  the 
absence  of  significant  sphere  wake  interference  effects  since  they  are  not 
included  in  the  theory.  At  higher  speeds  and  inclination  angles,  where 
sphere  wake  interference  would  be  expected  to  be  less  severe  and  not  a 
factor,  the  experimental  An  values  are  somewhat  greater  than  those  pre- 
dicted by  theory  and  the  B1B  values  are  in  very  good  agreement.  A possible 
explanation  for  this  may  be  in  the  less-than-perfect  Coleman  representation 
of  the  rotor  wake  longitudinal  "blow  back".  The  possibility  of  blade 
stall  tends  to  be  ruled  out  by  the  fact  that  the  average  blade  lift  coef- 
ficients are  approximately  the  same  value  across  the  lateral  axis  as  across 
the  longitudinal  axis. 

In  general,  the  predicted  cyclic  pitch  requirements  for  trim  are  in 
good  agreement  with  the  experimental  measurements,  and  it  can  be  concluded 
that  the  rotor  aerodynamic  representations  are  acceptably  accurate  for  the 
objectives  of  predicting  force  and  moment  equilibrium. 
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This  section  discusses  the  correlation  between  the  experimentally 
measured  transient  response  characteristics  of  the  AEROCRANE  model  in 

i 

various  configurations  and  flight  conditions  and  the  theoretical  pre-  j 

dictions.  Experimental  measurements  were  made  of  the  transient  response  j 

to  control  inputs  in  hovering  flight  at  two  center  of  gravity  positions  j 

with  various  sling  loads  and  levels  of  automatic  stabilization.  Table  II 
indicates  various  model  configurations  examined  in  hovering  flight. 

Figure  9 is  typical  of  the  data  obtained.  Additional  measurements 
are  presented  in  Figures  A-l  through  A-6  in  Appendix  A.  It  should  be  noted 
that  the  cyclic  pitch  trace  includes  the  cyclic  pitch  applied  by  the  auto- 
matic stabilization  system  as  well  as  that  applied  by  the  operator.  A j 

pulse  input  of  approximately  two  seconds  duration  with  an  amplitude  of  two 
degrees  was  employed  to  excite  the  transient  motion  of  the  model. 

The  transient  response  of  the  unstabilized  model  is  characterized  by  a 
lightly  damped  mode  with  a period  of  the  order  of  ten  seconds  and  a damping 

9 

ratio  of  less  than  0.1.  Figure  17  presents  the  frequency  and  damping 
characteristics  of  this  mode  as  measured  from  the  traces.  These  measurements 
should  be  viewed  as  approximations  to  the  character  of  this  mode  as  it  can  be 
seen  from  the  measured  time  histories  that  the  responses  are  not  precisely 
damped  sinusoids  indicating  the  presence  of  seme  of  the  other  modes  of  motion. 

Also  shown  is  the  comparison  between  the  theoretically  predicted  characteristics 
of  this  mode  as  given  in  Appendix  C with  the  measured  values.  The  agreement 
between  theory  and  experiment  is  excellent  for  the  high  center  of  gravity 
configuration.  The  damping  agrees  well  for  the  low  center  of  gravity  con- 

: 

figuration  and  the  frequency  agrees  for  run  11.  A lower  frequency  is  predicted 
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for  run  36.  The  reason  for  this  discrepancy  is  not  clear,  particularly  in 
view  of  the  excellent  agreement  for  the  other  three  runs. 

Figures  18  and  19  show  a direct  comparison  of  the  measured  and  calculated 
transient  response  characteristics  for  a cyclic  pitch  pul3e  input.  The  agree- 
ment in  general  is  very  good.  In  particular,  the  amplitude  of  the  roll  response 
predicted  agrees  very  well  with  the  measured  value.  The  theoretical  pitch 
response  exhibits  a smaller  amplitude  than  the  measured  response,  however, 
the  overall  agreement  is  very  good.  Thus,  the  theoretical  approach  presented 
here  gives  a good  prediction  of  the  vehicle  transfer  functions. 
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CONCLUSIONS 


Based  upon  the  analytical  and  experimental  research  program  reported 
herein  the  following  conclusions  are  made: 

1. )  The  Froude  scaled  model  of  the  AEROCRANE  vehicle  can 
be  flown  by  a remote  pilot  in  hovering  with  a resonable  level  of 
effort.  Remote  piloting  in  hovering  flight  was  made  considerably 
easier  by  incorporation  of  a relatively  simple  attitude  feedback 
system.  This  feedback  system  stabilizes  or  improves  the  damping 
of  the  unstable  or  lightly  damped  mode  characteristic  of  the 
AEROCRANE  in  hovering. 

2.  ) In  forward  flight  the  natural  damping  of  this  mode 
increases  and  no  difficulties  were  encountered  controlling  the 
vehicle  in  trimned  forward  flight  without  the  stabilization  system. 

3.  ) Accurate  quantitative  data  on  the  trim  and  dynamic 
stability  characteristics  of  the  model  were  obtained  frcm  remotely- 
piloted  flight  tests  in  a protected  environment. 

4. )  The  theoretical  model  for  the  forward  flight  trim  character- 
istics of  the  AEROCRANE  predicts  the  measured  experimental  data  with  a 
minor  adjustment  in  the  spherical  centerbody  drag  coefficient. 

5. )  The  analytical  model  of  hover  dynamic  stability  character- 
istics correlates  very  well  with  the  measured  model  transient  response 
indicating  that  the  analytical  model  provides  a good  representation 

of  the  transfer  functions  of  the  vehicle. 
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TABLE  I 


SUMMARY  OF  MODEL 

PHYSICAL  CHARACTERISTICS 

ROTOR  DIAMETER 

SPHERE  DIAMETER 

MODEL  WEIGHT  (WA) 

Low  c.g. 

39-9  ft. 

18.2  ft. 

177  lbs. 

High  c.g. 

19e  lbs. 

CENTER  OF  GRAVITY 

POSITION  (r  ) 

Low  c.g. 

2.71  ft. 

o 

High  c.g. 

1.83  ft. 

MOMENTS  OF  INERTIA 

1*  (pitch  and  roll)* 

Low  c.g. 

521  slug- ft 

RUNNING  WEIGHT  OF  UMBILICAL 

High  c.g. 

576  slug-ft 

653  slug-ft 

0.61  lb/ft 

CABLE 


* Determined  by  experiment  and  include  virtual  mass  effects. 
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SUMMARY  OF  HOVER  TEST  CONDITIONS 
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Ballast  is  rigidly  attached  to  gondola  base.  No  sling  load  is  used 


FORWARD  FLICfflT  TRIM  COHDITIONS 
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Figure  1.  Axis  Systems  Used  for  Analytical  Model. 

(a)  Gravity,  Body,  Rotor  and  Shaft  Systems. 
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Figure  1.  Continued. 

(c)  Load  and  Attacfciment  Point  Systems. 
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Figure  1.  Continued.. 

(d)  Model  and  Sling  Load  Free 
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Body  Diagram. 
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Figure  5.  Schematic  of  Attitude  Stability  Augmentation  System. 
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Figure  6.  Schematic  of  Angular  Rate  Stability  Augmentation 
System. 
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Theory:  Longitudinal  and  Lateral  Cyclic  Pitch. 
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Figure  17.  Period  and  Damping  of  Hovering  Retrograde  Oscillation, 
Comparison  of  Theory  and  Experiment. 
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APPENDIX  B 


UMBILICAL  CABLE  CONTRIBUTIONS  TO  DYNAMIC  STABILITY 

The  free -flight  model  of  the  AEROCRANE  was  equipped  with  an  umbilical 
cable  to  provide  power  and  control  signals  to  the  model  and  to  carry  data 
transducer  signals  to  the  ground  recording  instrumentation.  The  weight  of 
this  umbilical  cable  was  of  sufficient  magnitude  that  its  contributions  to 
the  dynamic  stability  of  the  vehicle  must  be  included  in  the  equations  of 
motion.  This  section  describes  the  manner  in  which  the  effects  of  the  um- 
bilical cable  were  incorporated  in  the  equations.  Only  the  hovering  flight 
case  is  considered  here. 

Figure  B-l  shows  the  configuration  of  the  umbilical  cable  for  hovering 
flight  experiments.  The  shape  of  the  umbilical  cable  is  assumed  to  be  a 
catenary.  The  lower  end  of  the  cable  is  in  contact  with  the  ground  and 
assumed  to  remain  fixed.  The  upper  end  of  the  cable  is  attached  to  the 
model  and  is  assumed  to  move  only  in  a horizontal  plane.  This  latter  assumption 

is  consistant  with  the  dynamic  stability  analysis  which  assumes  that  there  is 

e 

no  vertical  motion  of  the  model  during  its  transient  motion. 

The  cable  will  add  to  the  mass  of  the  vehicle  and  will  also  produce 
forces  and  moments  due  to  the  tension  in  the  catenary. 

The  effective  mass  of  the  cable  is  evaluated  first.  When  the  vehicle 
translates  it  will  not  carry  the  entire  mass  of  the  cable  with  it  but  rather 
will  deform  the  cable  into  a new  shape  as  shown  in  Figure  B-l.  It  is  assumed 
that  there  are  no  dynamics  associated  with  the  cable  itself.  That  is,  the 
motion  of  the  cable  is  assumed  to  be  quasi-static  such  that  its  shape  is 
always  that  of  a catenary.  The  effective  mass  of  the  cable  can  be  determined 
by  evaluating  the  kinetic  energy  of  the  cable  as  a function  of  the  horizontal 
translation  velocity  of  its  upper  end  (dx  ). 

O 
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Define  the  following  quantities 
m = running  mass  of  cable 


s « arc  length  of  cable 

x,y  = horizontal  and  vertical  coordinates  to  any  point  on  the  cable 
0(y)  = mode  shape  of  catenary,  normalized  by  horizontal  deflection  of 
upper  end  of  the  cable. 

Thus  the  equation  describing  the  local  horizontal  translation  of  the 


cable  is  given  by 


Ax  * Axq  0(y) 


(B-l) 


Differentiating  this  expression,  the  translational  velocity  of  the 
cable, incorporating  the  quasi-static  assumption  that  the  shape  of  the 
cable  is  unaffected  by  motion,  is  given  by 
Ax  * Axq  0(y) 

the  kinetic  energy  of  the  cable  is  therefore 


KE  » \ m J (Ax)3  ds 

Q, 

KE  * £ m Ax  8 J*  08  (y)  ds 


(B-2) 


Evaluation  of  the  integral  will  give  the  effective  mass  of  the  cable. 

Now  the  equation  for  the  perturbed  shape  of  the  catenary  is  developed 
such  that  the  integral  in  the  expression  for  the  kinetic  energy  can  be 
evaluated. 


The  equation  of  the  catenary  shape  shown  in  Figure  B-l  is  given  by 


y ■ LQ  (cosh  — - 1) 
o 


(B-3) 
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r; 

i 


whe: 


h = horizontal  component  of  tension  in  catenary. 

The  assumption  that  the  upper  end  of  the  cable  moves  only  in  a 
horizontal  direction  (Ay  - 0)  gives  a relationship  between  a pertur- 
bation of  the  tension,  Hq  (or  Lq)  and  a perturbation  of  the  upper  end 
of  the  cable,  Axq.  From  equation  (B-3) 

XXX  X 

Ay  = 0 = AL  (cosh  ' - l - -2  sinh  -2)  + Ax  sinh  [~]  = 0 (B-4) 

Lo  Lo  Lo  0 Lo 


Solving  equation  (B-3)  for  x 

x = L cosh  [1  + ] 

o L 

0 


(3-5) 


If  the  cable  translates  horizontally  a small  distance  Ax  equation  (B-5) 
becomes 

x + Ax  = (L  + AL)  cosh  [1  + - — x-77]  (B-6) 

O JLi  ■*"  All 

o 


t: 


Therefore  subtracting 'equation  (B-5)  from  (B-6) 


Ax  = (Lq  + AL)  cosh..  1 [1  + jr — y -^]  - Lq  cosh  1 [1  + ^-] 

o 


Assuming  that  the  change  in  tension  is  small  compared  to  the  initial  tension 


Ax  * AL 


cosh  (l  + •¥-)  - 
Lo 


JL 

"o 


Jv777! 


(B-7) 


Equation  (b-4)  is  used  to  eliminate  AL  from  equation  (B-7)  so  that  an 
equation  is  obtained  which  relates  the  local  horizontal  translation  of  the 
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i 


cable  to  the  translation  of  the  upper  end, 


sinhi  ^ 1 

4*  = 4x0  — -2 r jcoslf1  (1  + i) 

rslm  r + 1-C08tt  r 

o o o V 


f(f-)  * * (f) 

o o 


(B-8) 


Equation  (B-8)  thus  gives  the  mode  shape  of  the  perturbed  cable,  that  is, 
it  is  of  the  form  given  by  equation  (B-l)  such  that 


)sinh  ^ 

T F — | 

■=£■  sinh  — + 1 - cosh  — 
o o o 


cosh  1 (1  + 

o 


(B-9) 


Jv! + 2 (t’  | 

The  first  term  in  brackets  is  determined  by  the  initial  shape  of  the  cable, 
i.e.,  by  the  horizontal  distance  between  the  end  of  the  cable  on  the  floor 
and  the  end  of  the  cable  attached  to  the  model. 

It  is  possible  to  find  a simple  approximation  to  equation  (B-9)-  Figure 
B-2  shows  that  the  second  bracketed  term  in  equation  (B-9)  is  approximated 
very  closely  by 


k Jt 


* I cosh  (l  + t*£-)  - 


Kr1* 


(B-  10) 


Equation.  (B-9)  can  be  further  simplified  by  using  the  boundary  condition 
at  the  upper  end  of  the  cable  where  from  equation'  (B-l) 

0(^-1 

Lo 


Therefore 


XX  X 

~ sinh  ~ + 1 - cosh  -2. 
h L L 

0 0 0 


sinh  . — 


M 


(B-l!) 


Thus  the  approximate  mode  shape  is  given  by 


* <£>  -tk  ft 


(B-12) 


Substituting  equation  (B-l 2)  into  (B-2),  the  expression  for  the  kinetic 
energy  becomes 


KE  = i m Ax  * — J ds 
° yo  o Lo 


(B-13) 


The  arc  length  is  given  by  the  relationships  for  a catenary 


s = L sinh..  — 
o L 


(B-lU) 


Therefore  using  equations  (B-3)  and  (B-l4)  can  be  integrated.  The  result 
can  be  expressed  in  terms  of  arc  length  as,  noting  that  m sQ  = mc,  the  total 
mass  of  the  cable, 
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r 


In  (s  + ll  + s*)  

- 1 + 4 ( — M 


a 

o 


KE  = i m &&  * 
c o 


i 


1 + s. 


(B-15) 


where 


Lo,  the  ratio  of  the  horizontal  component  of  the  tension  to  the  running  mass 

of  the  cable  is  found  from  the  equation  of  the  catenary  (B-3)  knowing  the 

initial  end  points  of  the  cable.  For  the  various  hovering  experiments  the 

typical  height  of  the  model  (yQ)  was  80  feet  and  the  horizontal  distance  of 

the  model  from  the  point  at  which  the  cable  left  the  ground  was  25  feet  (x  ). 

o 

Substituting  these  values  into  equation  (B-3)  to  determine  Lq  and  then  using 
equation  (B-l4)  it  is  found  that, 
sq  = 10.79 

For  values  of  sq  of  this  magnitude  which  were  typical  of  all  of  the  hovering 
experiments,  equation  (B-15)  can  be  considerably  simplified.  For  sQ  greater 

f 

than  about  6 as  shown  in  Figure  B-2  a very  good  approximation  to  the  term 
in  brackets  in  equation  (B-15)  can  be  obtained  and  equation  (B-15)  can  be 
approximated  by 


KE  = im„ 


So 

2 

- 1 

i 

- 1 

o 

(B-l6) 


The  effective  mass  of  the  catenary  is  therefore 


m = m 

ce  c 


-2.  - i 


- 1 


(B-17) 


and  is  approximately  one-half  the  actual  mass  of  the  catenary.  The  value 
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■') 

i 

} 

one-half  was  used  in  the  dynamic  stability  analysis  since  the  cable  mass  I 

is  small  relative  to  the  mass  of  the  vehicle,  and  a more  refined  treatment  \ 

i 

was  not  considered  justified.  I 

\ 

Consequently,  the  inertia  effects  of  the  cable  are  considered  to  be  j 

represented  by  a concentrated  mass  equal  to  one-half  the  actual  mass  of 

the  cable  supported  (®ce)  located  at  the  attachment  point  of  the  cable. 

The  contributions  of  the  cable  to  the  acceleration  terms  in  the 

various  equations  of  motion  are 

Horizontal  Force 

AX  = m [u  + ZA  0] 
c ce  s 

Lateral  Force 

■ "ce  IV  2*  <B-l8> 

Pitching  Moment 

AMc  = ®ce  Z*  [V  Z*  ®] 

9 

Rolling  Moment 

ALc  = ttce  Z*  ['V  Z* 

The  terms  given  by  equations  (B-l8)  are  added  to  the  dynamic  stability 
equations  to  include  the  effect  of  cable  inertia. 

In  addition  to  these  acceleration  effects  the  weight  of  the  cable 
will  produce  forces  and  moments  on  the  vehicle. 


f , 


Figure  B-3  shows  the  geometry  involved  in  estimating  the  contributions 

of  the  cable  weight  to  the  equations  of  motion.  It  is  assumed  that  the 

cable  is  straight  when  looking  forward,  as  shown  in  the  figure.  The 

following  notation  is  employed , 

T =»  tension  in  cable  at  attachment  point 

=»  initial  slope  of  cable  at  attachment  point 

measured  with  respect  to  the  vertical 

AB  , dp  = perturbations  in  cable  slope  due  to  angular 
x y 

rotation  of  model. 

The  equilibrium  forces  applied  to  the  model  due  to  cable  tension  will 
result  in  an  initial  pitch  angle  (8i ) so  that 


- - T »*  * V 


T * 0 
c 

Z = T 
c 


(B-19) 


The  equilibrium  moments  due  to  cable  tension  are 


Mc  = - T ZA  Px 


L =»  0 
c 

The  perturbed  forces  and  moments  are 

Xc  + AX,  - - (T  + AT)(9i  + 0 + Bx  + A0x) 

Yc  + AYc  " (T  + + A®y) 

M + AM  * - (T  + AT)  z»  (9,  + 8 + + A*,) 


Lc  + ALc  = - (T  + AT)  Za  (0  + APy) 


(B-20) 
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The  slope  of  the  cable  can  be  found  by  differentiation  of  equation 


(B-3) 


J*  = sinh  ?- 
dx  L 

o 


(B-21) 


The  equation  of  the  catenary  gives  Lq  knowing  the  initial  position  of  the 


cable,  i.e.,  equation  (B-3)  is 


y = LQ  (cosh  ~ - 1) 


(B-3) 


Inserting  typical  values  from  the  hovering  experiments 
xQ  « 25  ft 
yo  = 80  ft 
Equation  (B-3)  gives 


L = 8.17  ft 
o 


consequently 

x_ 


r = 3.c7 
o 


= 9.79 

Lo 

For  these  typical  values  the  hyperbolic  functions  can  be  approximated  by 

x 

sinh  r~~  = ^ e 0 
o 

x 

.L 


cosh  7-  = t e 
Lo 


therefore  equations  (B-21)  and  (B-3)  can  be  approximated  as 


dz~  1 
dx  5 


x_ 

1 


(B-22) 


I 
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y » L0  (i  e 


(B-23) 


The  slope  given  by  equation  (b-22)  is  related  to  the  angle  Bx  by 


I - (90°- 


(B-24) 


Since  ?x  is  a small  angle 

dy  ~ 1 
to-9x 


From  equation  (B-23) 


(B-25) 


P = 2 e 
x 

Now  to  find  the  rate  of  change  of  Px  with  the  movement  of  the  end 
of  the  cable,  equation  (B-23)  can  be  expressed  in  terms  of  &x  as 


x - , x \ 2 

y ' (W > 11  r 


(B-26) 


Differentiating  equation  (B-26),  the  rate  of  change  of  Bx  with  x can  be 
found,  y is  constant  in  the  differentiation  since  it  is  assumed  that  the 
upper  end  of  the  cable  translates  horizontally  only.  The  result  is 


(i  - *g3 


y°l ln  rx  - (1  ’ ^ 


(B-27) 


Now  from  Figure  B-3  it  can  be  seen  that 
Ax  = Zk  AS 


Therefore  equation  (B-27)  can  be  written 


d0  y 


- »,>* 
f • a • 

x 


(B-28) 


Expressed  in  terms  of  the  initial  horizontal  position  of  the  cable 


<(,  !,  6* (1  - «*>  ^ «r 

49  ‘ *o  in  f - (1  - 6 ) 
px 


(B-29) 


Since  0x  is  assumed  to  be  small  this  result  may  be  further  approximated  by 


0x  lnr 

d»x  __  Z.  gx 

dF  " ^ ln  JL  - i 


(B-30) 


Using  typical  values 


xo  - 25  ft 


Z,  - 9 ft 


yQ  - do  ft 


t0  - 8.17 


Equation  (B-25)  gives  B 


B - 5.3 

x 

and  equation  (B-30)  gives 
<3B 

if  - o-01*? 

This  rate  of  change  of  the  cable  angle  with  vehicle  attitude  will  enter 
into  the  longitudinal  equations  of  motion. 


In  the  lateral  case  the  analysis  is  considerably  simpler  since  it  is 
assumed  that  the  cable  is  vertical  in  the  initial  condition.  Therefore, 
the  change  in  cable  angle  with  roll  is  a result  of  the  appearance  of  a 
component  of  the  initial  cable  angle  0x  in  the  lateral  plane.  Thus  from 
Figure  B-3 


<$„  Za 

__2l  = a 

d0  Px  xo 

For  the  typical  case  given  above 


(B-31) 


<*y 

rf  " °*°33 

The  variation  in  cable  tension  must  now  be  evaluated.  Since  the  cable 

tension  must  lie  along  the  direction  of  the  cable  and  the  vertical  component 

of  cable  tension  is  equal  to  the  weight  of  the  cable  supported, 

wc  ' 

T at  ^ 

cos  0x 

Therefore 


’x  cos3  0, 


sin  0 = W 0 

X C X 


(B-32) 


and  the  initial  tension  in  the  cable,  since  0x  is  a small  angle  is  equal 
to  the  weight  of  the  cable  T = Wc« 

Therefore  equations  (B-19)  can  be  written  as 


- wc  + 9i> 


Y - 0 
c 


(B-19) 


The  initial  cable  angle  will  produce  a small  inclination  of  the  vehicle 
in  hovering  equilibrium.  With  a buoyant  force,  ?a  , and  a sling  load 
weight  W* , the  complete  equation  for  x force  equilibrium  would  indicate 
an  initial  pitch  angle,  i.e.. 


(F,  - W*)  0±  - W^  (9i  +9X)  - 0 


9j 


Wc  ** 


F.  - Wtt  - V 


How  the  perturbation  forces  and  moments  can  be  evaluated  by  subtracting 
equations  (B-19)  from  equation  (B-20)  where 
' as 


AB, 


AB, 


dF 

AB 

_£ 

d« 


0 


T a W. 


dT  , 

41  - «rx  dF  8 


The  various  derivatives  are  given  by  equations  (B-30),  (B-31)  and 
(B-32). 


Therefore 

“c  - - "c  » + 9 ■ f.  ST  (,i  * >x>  9 


AB 

«c  - \ (1  * df>  * 


(8-33) 
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L ' 
I 

i ' 

I 


t 

I 


AMC  = - wc  z*  (i  + Jf)  e - Z,  ce±  + 3X)  0 

ALC  « - Z*  (1  - 5j)  0 


(B-33  Can't) 


dT 


Since  -jt-  ~ B and  -t-s-  ~ 3 the  second  term  in  the  horizontal  force  and 
dBx  x d 9 x 

pitching  moment  expressions  is  small  and  may  be  neglected. 

Thus  the  cable  weight  contributions  to  the  equations  of  motion  are 

«c  = - M=  ( 1 - df  > 9 


4I<=  = wc Cl  ♦ a i>  * 

4Mc  = - WoZ.  <1+  £f)  6 


(b-34) 


4Lc*-  V*  (1*r»)  » 

where  the  cable  angle  derivatives  are  given  by  equation  (B-30)  and  (B-31). 
These  terms  given  by  equations  (B-34)  taken  with  the  acceleration  terms 
given  by  equation  (B-l8)  thus  constitute  the  cable  contributions  to  the 
equations  of  motion. 

Since  for  the  typical  hovering  equilibrium  the  sample  calculation 

<*x  <*y 

presented  above  indicates  that  the  derivatives  and  are  small  com- 
pared to  1 an  average  value  equation  to  0.04  was  used  for  both  of  these 
derivatives  in  the  dynamic  stability  analysis. 

It  may  also  be  noted  that  the  cable  effects  described  in  this  section 
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give  rise  to  position  dependent  forces,  that  is,  a horizontal  translation 
of  the  vehicle  gives  rise  to  a translational  force  from  the  umbilical  cable. 

The  size  of  this  effect  may  be  estimated  by  calculating  the  trans- 
lational frequency  which  would  arise.  That  is  for  horizontal  translation 
the  equation  of  motion  is 


(m  + M ) Sc 


(B-35) 


where  it  has  been  assumed  that  the  cable  tension  is  equal  to  the  weight 
of  the  cable.  Thus  the  natural  frequency  arising  from  this  effect  is  from 
equation  (B-35) 


W 


c 

U'  + a^) 


d x 


The  the  typical  case  described  above 


0.0054 


For  the  low  center  of  gravity  experiments, 

t 

t 

m = 9.13  slugs 

and  a typical  cable  weight  lifted  is  50  lbs,  giving 
u>x  * 0.165  rad/sec 

giving  a period  of  38  seconds  which  is  considerably  slower  than  the  natural 
dynamics  of  the  vehicle  indicating  that  the  effects  are  small  and  may  be 
neglected. 
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APH2NDIX  C 


DETERMINATION  OF  AEROCRANE  MODEL  PHYSICAL  PARAMETERS 
AND  NUMERICAL  RESULTS  OF  DYNAMIC  STABILITY  ANALYSIS 

This  section  describes  the  manner  in  which  the  buoyant  force  and 

Inertial  properties  of  the  AEROCRANE  model  were  determined  experimentally 

and  also  presents  the  numerical  values  employed  in  the  stability  analysis 

as  well  as  the  results.  The  calculated  transient  response  characteristics 

based  on  the  numerical  values  given  here  are  discussed  elsewhere  in  this 

report. 

I.)  Determination  of  Buoyant  Force  and  Inertial  Properties 

Since  the  AEROCRANE  model  without  sling  load  and  umbilical  cable 
possessed  an  excess  of  buoyancy,  the  following  technique  was  used  to  deter- 
mine the  buoyant  force.  The  sling  load  was  placed  on  a scale  with  the 
model  floating  above  supporting  a length  of  umbilical  cable,  hQ.  The 
scale  reading  is  then  related  to  the  buoyant  force  by  the  following  equation 

[WQ  + WSL  + 0.61  hQ]  = S + Fb  (C-l) 

where 

WQ  = weight  of  the  basic  model,  lbs. 

WSL  = weight  of  sling  load,  lbs. 

S = scale  reading,  lbs. 

F*  * buoyant  force,  lbs. 

ho  = length  of  umbilical  cable  supported,  ft. 
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The  umbilical  cable  weighs  0.6l  lbs.  per  ft.  Table  C-l  lists  the  values 
of  the  buoyant  force  for  various  hovering  runs  determined  in  this  fashion. 
The  aerodynamic  thrust  is  determined  from  the  equilibrium  hovering  altitude 


r 


[ 


1*. 


' ? 


i 


for  the  specific  flight.  That  is, 

T - (WQ  + WSL  + 0.61  h^)  - F,  (C-2) 

where  T is  the  aerodynamic  thrust  and  hg^  is  the  measured  equilibrium 
altitude,  and  is  also  given  in  Table  C-l.  The  values  of  thrust  and 
buoyant  force  listed  in  Table  II  in  the  main  body  of  this  report  were 
determined  by  the  above  technique. 

The  moment  of  inertia  of  the  model  in  pitch  (equal  in  roll)  was 
determined  measuring  the  natural  frequency  of  the  model  oscillating  in 
pitch  with  the  rotor  not  rotating.  A weight  was  added  at  the  bottom  of 
the  model  to  increase  the  spacing  between  the  center  of  buoyancy  and  the 
center  of  gravity.  The  separation  of  the  center  of  buoyancy  and  the  center 

of  gravity  provides  a restoring  moment  and  consequently  produces  a free 

* 

motion  which  is  oscillatory  in  character.  The  restoring  moment  character- 
istic was  determined  experimentally  by  hanging  a weight  at  a blade  root 
and  measuring  the  angular  deflection  of  the  model  corresponding  to  this 
applied  moment.  The  equation  of  motion  for  the  pitch  oscillation  is 

I'  9 + Mg  9 * 0 (C-3) 

Consequently  the  natural  frequency  of  the  free  motion  is 
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The  measured  natural  frequency,  <%,  , and  the  restoring  moment  gradient. 

Mg , measured  as  described  above  were  used  to  determine  the  moment  of  inertia 
I1 . Note  that  this  procedure  determines  a moment  of  inertia  which  includes 
the  influence  of  accelerating  the  air  mass  adjacent  to  the  vehicle  (apparent 
mass  effects).  The  experimentally  determined  moment  of  inertia  in  pitch 
(equal  to  that  in  roll  because  of  symmetry)  is  listed  in  Table  C-l  as  well. 

To  obtain  the  configuration  of  runs  15  and  17  with  a higher  center  of  gravity 
position  a fifteen  pound  weight  was  added  at  the  upper  pole  of  the  spherical 
centerbody  and  this  is  reflected  in  an  increased  pitch  inertia. 

No  direct  measurement  of  the  polar  moment  of  inertia,  I , is  possible. 

z 

Consequently  this  value  was  calculated  based  on  knowing  the  size,  weight,  and 

location  of  all  of  the  various  components  in  the  model.  As  a verification  of 

this  procedure  the  moment  of  inertia  in  pitch,  i'  was  also  calculated  from 

component  contributions.  The  calculated  value  of  i'  was  within  one  percent 

of  the  value  determined  from  the  ocillation  tests  and  therefore  the  calculated 

value  of  I is  assumed, to  be  within  one  percent  of  the  actual  value, 
z 

Comparison  of  the  calculated  center  of  gravity  position  and  the  spacing 
between  the  center  of  gravity  and  center  of  buoyancy  which  can  be  determined 
from  the  measured  restoring  moment  characteristics  indicated  that  the  center 
of  buoyancy  was  0.39  ft-  above  the  plane  of  the  rotor  (the  geometrical  center 
of  the  centerbody  if  the  centerbody  is  a perfect  sphere)  indicating  a small 
distortion  of  the  centerbody  under  load.  This  small  distance  was  neglected 
in  the  analysis.  That  is,  the  center  of  buoyancy  was  assumed  to  lie  in  the 
plane  of  the  rotor. 
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Using  the  various  physical  parameters  of  the  model  given  in  Table  C-l 


and  the  equations  of  motion  presented  in  the  main  body  of  the  report,  the 
numerical  values  for  the  various  coefficients  in  the  equations  of  motion 
were  calculated  for  runs  11,  15,  17  and  36.  These  values  are  presented 
in  Figures  C-l  through  C-U.  It  should  be  noted  that  as  can  be  seen  by 
comparison  of  the  numerical  values  of  the  elements  of  the  damping  and  control 
matrices  with  the  literal  expressions  given  on  pages  35  and  37  that  the  rotor 
inplane  forces  were  neglected.  Previous  analyses  have  shewn  that  these 
terms  have  only  a minor  effect  on  the  dynamic  stability  and  response 
characteristics  in  hover. 

Table  C-2  presents  the  eigenvalues  calculated  for  these  four  hovering 
cases.  Also  presented  are  the  corresponding  period  and  damping  ratio  as 
well  as  an  identification  of  each  mode.  The  dynamic  motion  of  the  AEROCRANE 
model  in  hovering  in  the  configuration  flown  is  characterized  by  five 
oscillatory  modes.  Ttjree  of  these  modes  are  associated  with  the  vehicle 
and  two  with  the  sling  load.  One  of  the  basic  modes  of  the  vehicle  is  a 
fast  motion  that  is  well  damped  with  a period  of  the  order  of  1.8  seconds 
and  a damping  ratio  of  the  order  of  0.5.  This  is  basically  an  angular 
motion  with  its  character  determined  primarily  by  the  aerodynamic  damping 
of  the  rotor  and  the  gyroscopic  moments.  There  is  a lightly  damped  mode 
t with  a period  of  the  order  of  11  seconds  and  a damping  ratio  of  0.1  or 

less.  This  is  the  mode  which  dominates  the  measured  transient  response 
characteristics  owing  to  its  small  damping  ratio.  As  discussed  in  Ref- 
erence 1,  this  mode  may  be  described  as  a retrograde  mode,  that  is,  as 
! a result  of  the  polar  symmetry  of  the  vehicle  this  transient  motion  is  in 
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fact  a circling  motion  of  the  vehicle  and  in  this  mode  the  circling  takes 
place  opposite  to  the  direction  of  rotor  rotation.  The  remaining  vehicle 
mode  has  a period  of  the  order  of  20  to  30  seconds  and  is  well  damped 
with  a damping  ratio  of  the  order  of  0.7.  It  can  be  characterized  as 
an  advancing  mode  as  it  corresponds  to  circling  in  the  direction  of  rotor 
rotation.  Owing  to  its  large  damping  ratio  the  presence  of  this  mode  is 
not  apparent  in  the  measured  or  calculated  transient. responses.  The 
remaining  two  modes  are  associated  with  the  sling  load  motion  in  two 
directions.  The  period  of  these  motion  is  of  the  order  of  4 seconds 
and  the  damping  ratios  are  very  small  as  the  sling  load  damping  was 
neglected.  The  isolated  sling  load  period  is  4.69  seconds  indicating 
that  there  is  sane  coupling  with  the  vehicle  motion. 

Comparison  of  the  calculated  transient  response  based  on  these 
numerical  values  with  the  experimentally  measured  transient  response 
is  discussed  in  the  main  body  of  the  report. 
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EIGENVALUES 

sec”3" 

PERIOD 

sec. 

DAMPING 

RATIO 

MODE 

-2.275  ± 3.575i 

1.76 

0.536 

FAST  ANGULAR  RESPONSE 

-0.086  ± 0.668i 

9.41 

0.128 

LIGHTLY  DAMPED  RETROGRADE 

-0.248  ± 0.2721 

23.10 

0.674 

WELL  DAMPED  ADVANCING 

-0.035  ± 1.6151 

3.89 

0.022 

SLING  LOAD 

-0.093  ± 1.4381 

4.37 

0.065 

SLING  LOAD 

-2.096  ± 3.3711 

1.86 

FAST  ANGULAR  RESPONSE 

-0.013  ± 0.4841 

12.98 

LIGHTLY  DAMPED  RETROGRADE 

-0.208  ± 0.1931 

32.55 

0.733 

WELL  DAMPED  ADVANCING 

-0.004  ± 1.4111 

4.45 

0.003 

SLING  LOAD 

-0.ce9  ± 1.3651 

4.60 

0.021 

SLING  LOAD 

RUN  NO.  17. 

9 

-2.300  ± 3-3441 

1.88 

0.567 

FAST  ANGULAR  RESPONSE 

-0.008  ± 0.5931 

10.60 

0.013 

LIGHTLY  DAMPED  RETROGRADE 

-0.244  ± 0.3311 
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-0.013  ± 1.5041 

4.18 

0.009 

SLING  LOAD 

-0.069  ± 1.4001 

4.49 

0.049 

SLING  LOAD 

Ln  no.  36. 

-2.361  ± 3.8341 

1.64 

0.524 

FAST  ANGULAR  RESPONSE 

-0.051  ± 0.5641 

n.i4 

0.090 

LIGHTLY  DAMPED  RETROGRADE 

-0.229  ± 0.2101 

29.90 

0.737 

WELL  DAMFED  ADVANCING 

-0.015  ± 1.5381 

4.08 

0.010 

SLING  LOAD 

-0.069  ± 1.4151 

4.44 

0.049 

SLING  LOAD 

t I 
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Figure  C-l:  Numerical  Values  of  Matrix  Elements  for  Run  No.  11. 
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Figure  C-2:  Numerical  Values  of  Matrix  Elements  for  Run  No.  17. 


117 


MASS  MATRIX  (M) 


1.089 

0 

0 

0 

0 

0.0068 

0 

I.O89 

0 

0 

-0.0088 

0 

1.5 

0 

1 

0 

0 

-0.0555 

0 

1.5 

0 

1 

0.0555 

0 

0 

-0.5108 

0 

0 

1.062 

0 

0.5108 

0 

0 

0 

0 

1.062 

DAMPING  MATRIX  (O 

3.03 

3.66 

0 

0 

-0.0379 

0.1543 

-3.66 

3.03 

0 

0 

-0.1543 

-0.0379 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

SPRING  MATRIX  (k) 

1.709 

0 

-0.881 

0 

0 

0 

0 

1.709 

0 

-0.881 

0 

0 

1.79 

0 

1.79 

0 

0 

0 

0 

1.79 

0 

1.79 

0 

0 

0 

1.4 

0 

-5.59 

0 

0 

-1.4 

0 

5.59 

0 

0 

0 

CONTROL  MATRIX  (F) 

0 

-8,82 

8.82 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Figure  C-3:  Numerical  Values  of  Matrix  Elements  for  Hun  No.  36. 
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Figure  C-4:  Numerical  Values  of  Matrix  Elements  for  Run  No.  15 


The  equations  for  equilibrium  flight  formulated  in  the  text  of  the 
report  do  not  include  the  sling  load  or  the  umbilical  cable.  These  effects 
can  be  incorporated  in  the  following  fashion. 

Since  the  aerodynamic  forces  on  the  sling  load  are  negligible  owing 
to  the  high  density  of  the  load  and  the  low  velocity  of  the  experiments 
the  equilibrum  flight  orientation  of  the  sling  load  cable  is  vertical  as 
can  be  seen  from  the  steady-state  solutions  of  the  latter  two  of  equations 
(65).  Thus 

0L  * -6 
0L  a -0 

Consequently,  as  seen  from  the  first  four  of  equations  (65)  the  sling  load 
adds  to  the  weight  in  the  two  force  equations  and  produces  a moment 
about  the  center  of  gravity  of  the  vehicle.  If  it  is  assumed  that  the 
angle  that  the  umbilifcal  cable  makes  with  the  vertical  is  negligible  in 
forward  flight,  the  umbilical  also  appears  in  the  equilibrium  equations  in 
a similar  fashion  to  the  sling  load,  as  may  be  seen  from  equations  (67^. 

The  contributions  of  the  umbilical  from  equation  (67)  and  the  sling  load, 
from  equation  (65)  to  the  force  balance  equations  given  by  equation  (8) 

• are 

horizontal  force 

(-mtg  - WQ)  sin  9 cos  0 
side  force 

(mtg  + WQ)  sin  0 
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vertical  force 

(mug  + Wc)  cos  0 cos  S 

Therefore  for  the  force  balance  equation^  to  include  the  effects  of  the 
sling  load  and  the  umbilical  cable,  the  term  W may  be  interpreted  as 
W = WQ  + mtg  + Wc 

That  is,  it  is  the  total  weight  supported  by  the  buoyant  force  and  aerodynamic 
thrust.  This  quantity  W is  referred  to  as  the  gross  weight  elsewhere  in  the 
report. 

For  inclusion  in  the  mcoent  equilibrium  both  of  these  weights  act  at 
their  attachment  point  a distance  ZA  below  the  center  of  gravity  of  the 
vehicle. 

Thus  to  equations  (13),  the  following  terns  must  be  added 
pitch  moment 

ZA  C“l8  + sin  0 cos  0 

c 

roll  moment 

♦ 

-ZA  [mtg  + Wfi]  sin  0 

Thus  to  sunmarize,  to  include  the  effects  of  the  umbilical  cable  and  the 
sling  load  in  the  trim  calculations,  the  term  W can  be  interpreted  in  the 
force  balance  equations  (equation  8)  as  including  the  weight  of  the  sling 
load  and  the  umbilical  cable.  The  moment  equilibrium  equations  (equation 
16)  must  be  modified  to  include  the  above  terms.  The  simpler  form  given 
by  equation  (18)  can  not  be  used. 
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